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The  Workshop  on  Nonlinear  Interactions  in  Magnetic  and  Magneto¬ 
optic  Materials  was  held  December  12-14, 1993,  at  the  Cotintryside  Inn  in 
Costa  Mesa,  California,  It  was  sponsored  by: 

California  Coordinating  Committee  for  Nonlinear  Studies 
Institute  for  Surface  and  Interface  Science,  University  of  California, 
Irvine 

Physics  Department,  University  of  California,  Irvine 
U.S.  Army  Research  Office. 


The  organizing  committee  was: 

R.  F.  Walhs  (Chair),  D.  L.  Mills,  C.  S.  Tsai;  U.  C.  Irvine 

R.  F.  Soohoo;  U.  C.  Davis 

A.  D.  Boardman;  University  of  Salford. 

The  Workshop  was  attended  by  42  participants  representing  the 
United  States,  United  Kingdom,  Russia,  China,  Italy,  Japan,  Germany,  Israel 
and  Brazil.  The  program  consisted  of  21  invited  talks  and  6  contributed 
talks.  Topics  covered  were  critical  behavior  and  phase  traditions  in 
magnetic  systems,  nonlinear  spin  waves  and  magnetic  solitons,  properties 
and  device  applications  of  magnetooptic  materials,  magnetic  resonance  and 
chaos  in  nonlinear  magnetic  materials,  and  giant  magnetoresistance.  The 
detailed  program  follows: 

Saturday  12/11/93 

6:00  -  7:30  PM  Registration 


Sunday  12/12/93 

8KX)  -  8:45  AM  Registration 

8:45  •  9:00  Welcome  address  and  administrative  information 

Session  1  Chair:  P.  E.  Wigen 


9:00-9:40  H.  Suhl,  UC  San  Diego 

"Critical  Behavior  of  Magnetic  Domains 


9:40  - 10*^0 

10:20  - 10:50 
10:50  - 11:30 


11:30  - 12:10 


D.  L.  MiUs,  UC  Irvine 

"Field  Induced  Phase  Transitions  in  Magnetic  Superlatnces 
Coffee 

S.  Nikitov,  Institute  of  Radioeng.  &  Elect.,  Moscow 
"Nonlinear  CW  Magnetostatic  Spin  Waves  in  Ferromagnehc 
Films:  Modulation  Instability,  Decay  Processes,  Multistabihty, 
and  New  Effects" 

A- D.  Boardman,  U.  of  Salford  ,  „  ^ 

"Fundamental  Theory  of  Bright  and  Dark  Solitons  in  Magnetic 

Media" 
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12:10  -  2:00 

Lunch 

Session  2  Chair:  W.  J.  Miceli 

2:00-2:40 

C.  S.  Tsai,  UC  Irvine 

"Magnetostatic  Waves-Based  Guided  Wave  Magnetooptic 
Interactions,  Devices  and  Applications" 

2:40  -  3:20 

F.  HeUman,  UC  San  Diego 

"Anisotropy  in  Amorphous  Materials  for  Magnetooptics" 

3:20  -  3:50 

Coffee 

3:50  -  4:30 

H.  H.  He,  Huazong  Univ.  of  Science  and  Tech.,  China 

"Magnetooptic  Thin  Films  and  Waveguide  Isolators" 

4:30-5:10 

Y.  Miyazaki,  Toyohashi  Univ.  of  Technology,  Japan 
"Switching  and  Filtering  Characteristics  of  Optical  Devices 
Using  MSSW  in  Multilayered  Magnetooptic  Waveguides" 

5:10  -  5:50 

Y.  Fetisov,  Inst,  of  Radioeng.,  Electr.  and  Autom.,  Moscow 
"Anomalous  Optical  Waveguide  Modes  Conversion  in  YIG  Film 
under  Microwave  Pumping" 

Monday  12/13/93 


Session  3  Chair:  C.  E.  Patton 

9:00  -  9:40 

P.  E.  Wigen,  Ohio  State  U. 

"Controlling  Chaos  in  Nonlinear  Magnetic  Thin  Films" 

9:40  - 10:20 

S.  M.  Rezende  and  F.  M.  de  Aguiar,  Fed.  U.  of 

Pernambuco,  Brazil 

"Crises  and  Control  of  Chaos  in  Spin  Wave  Instabilities" 

10*.20- 10:50 

Coffee 

10:50- 11*.30 

R.  F.  Soohoo,  UC  Davis 

"Internal  Fields  in  the  Ferromagnetic  Resonance  of  Soft 
Ferromagnetic  Materials" 

11:30-12:10 

A.  Slavin,  Oakland  University 

"Sample-Size  Effect  in  Spin  Wave  Auto-oscillations" 

12:10  -  2:00 

Lunch 

Session  4  Chair:  M.  Ciftan 

2:00  -  2:40 

P.  Kabos,  Colorado  State  U. 

"Paramagnetic  Magnon  Wavevector  Distribution  in  Yttrium  Iron 
Garnet  Films  by  Brillouin  Light  Scattering" 

3 


2:40  -  3:20 

R.  Marcelli  &  P.  De  Gasperis,  Inst,  of  Electronics,  Rome 
Tower  Thresholds  and  Frequency  Shifts  in  the  Excitation  of 
Magnetostatic  Volume  Wave  Solitons  in  Delay  Line 
Configurations'* 

3:20  >  3:50 

Coffee 

3:50  -  4:30 

B.  A.  Kalinikos,  St  Petersburg  Electrotechnical  Univ. 

"Spin  Wave  Envelope  Solitons  in  YIG  Films:  Formation, 
Propagation,  Reflection  and  Collision  Effects" 

4:30  •  5:10 

D.  Young,  Rockwell  IntL  and  C.  S.  Tsai,  UC  Irvine 

"Some  Observations  of  Nonlinear  Phenomena  in  Magneto¬ 
static  Waves  and  Guided-Wave  Magnetooptic  Experiments" 

5:10  >  5:30 

K.  Booth,  U.  of  Salford 

"Observation  of  Spatial  Solitons  in  Magnetic  Thin  Films 

5:30  >  5:50 

M.  Chen,  Colorado  State  U. 

Backward  Volume  Wave  Microwave  Magnetic  Envelope  Solitons 
in  Yttrium  Iron  Garnet  Films" 

7K)0 

Banquet 

Tuesday  12/14/93 

Session  5  Chair:  R.  Marcelli 

9:00-9:40 

A.  Berkowitz,  UC  San  Diego 

"Giant  Magnetoresistance  and  Intercluster  Exchange 
Interactions  in  Heterogeneous  Alloys" 

9:40  •  10:20 

R.  E.  Camley,  U.  of  Colorado,  Colorado  Springs 

"Giant  Magnetoresistance  in  Magnetic  Multilayers" 

10:20  - 10:50 

Coffee 

10:50  - 11:30 

I.  Schuller,  UC  San  Diego 

"Connection  Between  Structure  and  Giant  Magnetoresistance" 

11:30  - 12:10 

A.  C.  Ehrlich,  Naval  Research  Laboratory 
"Quantum-Well  States  and  Giant  Magnetoresistance  in 
Magnetic  Multilayers" 

12:10  -  2:00 

Lunch 

Session  6  Chair:  A.  Slavin 

2:00  -  2:20 

M.  Tsankov,  Colorado  State  U. 

"Forward  Volume  Wave  Microwave  Magnetic  Envelope  Solitons 
in  Yttrium  Iron  Garnet  Films 
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2:20  ~  2:40 


2:40  -  3:00 


3:00  -  3:20 


3:20  -  3:50 
3:50  -  4:50 


4:50  -  5:00 


A.  Pomyalov,  Bar  Ban  University,  Israel 

"The  Fine  Structure  of  Subsidiary  Absorption  at  Parametric 
Excitations  of  SW  in  YIG  Films" 

M.  M.  Shabat,  Islamic  University  of  Gaza 
"Nonlinear  Electromagnetic  Surface  Waves  Guided  by  a 

Groimded  Gyromagnetic  Film" 

N.  Waby,  U.  of  Salford 

"On  the  Existence  of  Envelope  Sobtons  in  Antiferromagnetic 
Thin  Films" 

Coffee 

Panel  Discussion  on  Imjiortant  Issues  and  Future 
Prospects 

Moderator:  D.  L.  Mills 

Panelists:  A.  Berkowitz,  A.  D.  Boardman,  C.  E.  Patton, 
C.  S.  Tsai,  P.  E.  Wigen 

Closing  Remarks 


The  Workshop  served  a  very  worthwhile  function  in  bringing  together 
experts  from  all  over  the  world  in  a  setting  conducive  to  intense  one-on-one 
interactions.  Particularly  significant  was  the  interaction  between 
experimentalists  and  theorists  which  led  to  a  clarification  of  important  issues 
concerning  magnetic  solitons.  Also  of  great  interest  was  the  overview  of 
present  and  prospective  applications  of  giant  magnetoresistance  and 
magnetooptic  materials  to  useful  devices. 
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■  Nonlinear  effects  are  important  in 
science  and  engineering  , 


■  They  can  be  used  to  transmit  and 
control  pulses  in  various 
waveguides 


■  Frequency  sweep  across  pulse 
envelopes 

■  Phase,  amplitude,  frequencies  are  • 
changed 
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■  fisherman’s  warning 

■  seventh  wave  wiii  sink  you 


SECH 


ENVELOPE 

7th  WAVE  WILL  SINK  YOU 
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The  Nonlinear  Schrodinger  Equation 

•  ■ _ 

r  '  '  .  ■  •  '  , 

H  Maxwell  +  Landau-Lifshitz  equations 

■  '  .  « 

•  .  equations  for  amplitudes  of  magnetic  potentials  q/ 

•  H=-Vvj/ 

^  application  of  boundary  conditions 

D(w,kjAh.A.=  0 
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Nonlinear  Wave 
e^ansion  about  (k 


The  Nonlinear  Schrodinger  Equation 

■  A  is  the  complex  amplitude  of  the  magnetic  potential 

I  ' 

■  Angular  frequency  (©),  wave  vector  (k) 

■  Choose  a  working  point 

(O)o,ko) 

■  weakly  nonlinear  weakly  dispersive 

^  expand  =  0  about  ((Oo,ko,0) 

■  Define  small  deviation  K  from  ko 

K=k-ko 

i  .  ,  •  • 

ko  =  (0,0.ko) 

K=  (Kp  ,  0  ,  Kz) 


II 


Wavenumber 


I 


.  dA  P2  I  ^|2  /  ^ 


f 


T  =  t-zP, 


/■  .  A 

V/ 


T  measured  in  frame  of  reference  moving  with  pulse 
d^k  1  3^00 


P2  = 


5® 


dk 

Y  i  .  I  ,2 


1  9® 


Q\A\ 


rp2 

L  — 


:  dispersion  length 
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:  nonlinear*  length 


no 


DIMENSIONLESS  FORM  TEMPORAL 


.dA  0)2  d'^A 
2 


0)2  = 


Defme^  =  7f'/ 
Hence  [s'=§  v 


‘^2  . 
^  =  -::rs 


i^^^^^^sgnirJP\Dl\4A 


A=  ^u- 

define  'y  g  ’  Po  is  incident  power 

Lighthill  criterion 
drop  thp  primes  4* 


,-^,i^+it/rt/=o 

2  <7^  '  ' 


^2_ 


2  02S  U-Nu 
N  is  the  soliton  order 


=  0^ - 

“  FIBRE 

uCzi.tjV,)* 

U(zi+i  ^  tj+i ) 

/^^U(Zi*.  tj) 

U(zi,tj)*  • 

U(Zi,tH  )• 

U(Zi+,  tj.i) 

• 

*■  z  =  L 


■  we  can  use  u(zi  tj+i)  uCz,  t,., )  and  U(zi,tj) 

to  find  the  next  step  along  the  fibre  U(zi+i  tj) 


JZI 


Modems-Telephone 
Service  for  Computers 


a.  Carrier 


0  01101  00010 


b.  Amplitude  Modulation 

0011010  0010 


c.  Frequency  Modulation 


d.  Phase  Modulation 


r 


NO  SOLITARY  NAVES  ^  SOUTARY  lAVES 


•  define  a  aurface  D^-f(k) 


3rd- Order  Soliton 


X'|Tsua'|ui 


0.000  -4.0 


Soliton  Collision 


opn'|T|duiv 


axis  I  axis 


axis 


Distance 


0»li-r  i  Wi,  +  1  u  1  =  0 

i  is  'propagation'  distance 

\  'tfi  is  transverse  variable 


;r>i^nf|arY  rnflditlonS: 


,2 

U  I  « 


•tir 


(xV 


along  whole 

guided^  f  . .  . 


'  .2  ♦ 
H . / 


causal  violation 


. : . . vm  ■ . 

can  use  beam  agation 


SPIN  WAVE  ELECTRONICS 


80LIT0N  FORMATION 
DECAY  INSTABILITIES 

SELF-INTERACTION 
MAGNETOELASTIC  WAVES 

BISTABILITY 

CONSTRUCTION  OF  rf  DEVICES 
PARAMETRIC  INTERACTIONS 


MSW  Delay  Line 


ABOUT  DC  FIELD 


(COh+COhCOm) 


Zvezdin,  Popkov  (1983) 


Forward 


- -k  coH=T(Hi4irM,)  (0,=4tTM^ 

Power  Increase  Mzi  cOhT 
Therefore  positive  nonlinear  change 


6  k' 


<0 


A  nl  - 

Ao)  >0 


(co^+(o„co„f  Backward 


00M=4irTM 


k  k  2 

M.i  (oa  ^,>0  aJ'<0 

z  ok 


COh+^J^m 
— .  2 


Surface 


,,  2  ,  , ,  , ,  Y////////ZA^//Z/Z/V-//A 

((Oh+Wh^m)  2 

M  i  a)..i  ^,<0  A(o‘''<0 


Data  used: 


,  ■  width  ofsquare  input  pulse  0.18/1^ 

j 

6'CJOOl 

■  normalised  input  intensity  =  0.0009 

; .  0*0< 

I 

\ 

;  ■  group  velocity  of  linear  wave  packets  ~  1-8  cm  / 

1 ' 

‘  ■  dispersion  of  medium  co^  =  -0.3  cmV \is 

t 

I 

■  derivative  of  the  spin  wave  frequency  with  respect 
. to  \Uf  =30,000/|LLy 


■  damping  coefficent  'H  =  6.0 


43 
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Damped  Pulse  On  A, 5  Micron  YIG  Film 
Square  Pulse;  Inlilit  Amplitude  0.1 


axis 


Undamped  Pulse  On  A  Thin  YIG  Film 
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The  Critical  Condition 
of  Self-trapping  (Self-focusing) 


4 :  ray*  ratum  to  axia 


S'3 


focusing 


Propagation  in  a  Waveguide 


U  take  advantage  of  waveguide, 
power  at  balancing  condition: 


xrd 


8a  r 


d  is  guide  dimension,  focusing  changes  r  but  not  d 


P<R 


P>B. 


decreases 


cross  P=lg 


over 


cross  P=lg 
over 


l_  ■  '*[  *  ^  •  i|  I  I  *  •  .  '  ^  »  ; 


.1  ,  ‘ 


4  P=R  ,  balance  obtained 


S'S 


otSiHIo  (f%f%  P'8 » 


Propagation  in  a  Waveguide 


take  advantage  of  waveguide, 


d  IS  guide  dimension,  focusing  changes  r  but  not  d 


P<R. 


P>P« 


diffraction 
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A  model  which  can  account  for  the  experimentally  observed  variations  of  the  giant  magnetoresistance 
in  thin  magnetic  multilayers  with  mean  free  path,  interface  roughness,  magnetic  layer,  and  normal  layer 
thickness  has  been  developed.  The  model  requires  the  existence  of  quantum-well  states  within  individual 
layers  or  groups  of  layers,  depending  on  the  magnetic  state  of  the  film.  The  calculated  results  are  ob¬ 
tained  by  the  application  of  quantum  size  effect  transport  theory  to  these  indivtdual  layers. 


PACS  number?:  73.50JI,  73. 61. At,  75.70.H 

The  first  report  by  Baibich  et  j/.  [ll  of  a  very  large 
negative  magnctorcsistancc.  the  so-called  giant  magne- 
toresistance  (OMR),  in  Fe/Cr  multilayers  has  provoked  a 
great  deal  of  experimental  work  in  sandwiches  and  multi¬ 
layers  of  this  system  l2.3l  as  well  as  other  magnetic  mul¬ 
tilayers  t4l.  This  large  decrease  in  resistance  with  ap¬ 
plied  magnetic  field  occurs  when  an  originally  antiparallcl 
orientation  of  the  magnetization  of  adjacent  ferromagnet¬ 
ic  layers  of  the  film  is  driven  parallel.  The  origin  of  the 
magnetic  structure  and  behavior  of  these  systems  is  closc- 
*  ly  related  to  that  of  the  GMR  and  like  it  is  also  under  in¬ 
tense  study.  Currently,  experiments  are  exploring  and 
finding  dependence  of  the  properties  of  these  multilayers 
on  the  detailed  electronic  structure  of  the  constituent  ma¬ 
terials  15).  In  this  paper  only  the  behavior  of  the  magne¬ 
toresistance  arising  from  the  relative  orientation  of  the 
magnetization  will  be  addressed.  We  will  use  AC  and 
“3V”  to  mean  the  magnetic  and  nonmagnetic  layers,  re¬ 
spectively.  of  any  of  these  films. 

The  goal  of  this  paper  is  to  define  a  simple  and  plausi¬ 
ble  mechanism  which  can  semiquantitatively  account  for 
the  large  number  of  well-established  and  still  puzzling 
features  associated  with  these  multilayered  films.  These 
include  the  following;  a  larger  effect  in  a  superlatiice  than 
in  a  sandwich  structure;  a  larger  effect  at  low  tempera¬ 
tures  than  at  room  temperature;  a  clear  enhancement  of 
the  GMR  with  interface  roughness  l6l;  a  rapid  variation 
of  the  effect  with  N  thickness  l7l;  a  slower  variation  with 
A/  thickness  which  can  have  either  of  two  characteristics 
(4l,  a  monotonic  decrease  of  the  GMR  with  Af  layer 
thickness  above  some  relatively  small  thickness  (10  A  in 
Fe/Cr)  or  an  increase  of  GMR  up  to  large  thicknesses  of 
the  magnetic  layers  followed  by  a  broad  maximum  as  the 
thickness  is  increased  further,  and  a  behavior  dependence 
on  the  detailed  electronic  structure  of  the  multilayer  ma¬ 
terials.  The  model  we  present  here  is  a  single  mechanism 
which  predicts  the  existence  of  all  of  these  features  in  a 
very  natural  way. 

The  theoretical  modeling  has  generally  emphasized 
spin  dependent  scattering  l2.4l.  More  specifically.  Inoue. 
Oguri,  and  Maekawa  I8l  addressed  the  origin  of  the  spin 
dependent  potential  responsible  for  the  spin  dependent 
scattering  at  the  M/N  interface.  Camley  and  Barnas  l91 


and,  with  fewer  approximations.  Hood  and  Falicov  [lOl 
have  extended  Fuchs'  semiclassical  theory  (ill  of  size 
effect  to  spin  dependent  scattering  at  the  Fe/Cr  interface. 
Levy.  Zhang,  and  Fert  (121  have  followed  the  quantum 
mechanical  theory  of  Tesanovic.  Jaric.  and  Maekawa 
(13)  to  treat  the  Fe/Cr  interface  as  a  source  of  spin 
dependent  intcrfacial  roughness  scattering.  Stearns  [14] 
has  treated  the  spin  dependent  scattering  as  arising  from 
the  difference  in  the  spin  up  and  spin  down  density  of 
states. 

In  this  paper  it  is  shown  that  the  GMR  can  be  largely 
understood  in  terms  of  (i)  the  vacuum-metal-like  charac¬ 
ter  of  the  interface  between  A/  and  N  and  the  consequent 
quantum-well-like  states  established  within  each  layer  or 
groups  of  layers  of  certain  multilayer  films  and  (ii)  the 
application  of  already  existing  quantum  mechanical  size 
effect  theories  of  electrical  conduction,  (i)  has  already 
been  experimentally  established  in  a  number  of  different 
overlayers  of  nonmagnetic  metals  on  other  nonmagnetic 
metal  substrates.  Furthermore.  Edwards  and  Maihon 
(15.161  have  presented  a  theory  of  the  exchange  coupling 
of  the  ferromagnetic  layen  in  magnetic  multilayers  based 
on  the  idea  of  quantum  confinement  of  the  electron' 

within  the  individual  layers. 

Recently,  using  angle-resolved  photocmission  Lindgrc 
and  Wallden  (17)  and  Chiang  and  co-workers  (18.19 
have  discovered  and  explored  discrete  valence  electro 
states  in  thin  metal  overlayers  on  a  metal  substrate 
These  have  been  identified  as  arising  from  electro 
confinement  within  the  overlayer  and  thus  as  quantum 
well  states  (QWS)  (18).  The  degree  to  which  the  ovei 
layer  states  manifest  quantum-well  behavior  revolve 
around  the  amount  of  electron  reflection  versus  electre 
transmission  ll8l.  This  depends  on  the  dissimilarity  . 
the  two  metals  at  the  interface  and  the  abruptness  of  tr 
interface.  Nevertheless.  Ag(l  1  D  on  Au(lll).  two  ve 
similar  metals,  displays  strong  QWS  characteristics 
those  Ag  electrons  energetically  below  the  Ag  L  point  ar 
above  the  Au  L  point  in  the  Brillouin  zone  boundary, 
this  region,  called  by  Chiang  and  co-workers  (18)  a  ••re! 
live”  gap.  there  arc  no  energy  gaps  in  either  material,  b 
there  are  no  available  states  in  Au  which  conserve  bo 
energy  and  the  (111)  component  of  the  Ag  electror 
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Bases  of  the  Model 


The  interfaces  in  these  films  can  give  rise  to 
quantum  well  states, 

[Chaing  et  al,  Lindgren  and  Wallden 
(photoemmision),  Himpsel(inverse  photoemission)] 


The  band  structure  of  the  non-magnetic 
material  is  much  more  like  one  of  the  spin 
bands  of  the  magnetic  material  than  the  other. 
(Inoue  et  ai.,  Papaconstantopoulos  (Fe-Cr)) 


The  dimensions  of  the  wells  depend  on  ferro  or 
antiferromagnetic  orientation  of  the  magnetic 
layers. 


Thus  a  magnetoresistance  arises  from 
established  quantum  size  effect  considerations. 
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EXTENSION  OF  MODEL  TO: 


unequal  mean  free  paths  in  the  different  layers. 

a  much  smaller  electron/atom  number  (0.131). 

various  numbers  of  layers  (sandwiches  to 
superlatices  of  63  pairs  +1). 


r_L  L  L  L  L  t  L  t  t  I  t  L  t  tit  t  L  i_  L  L  I  I  I  111  t  I  L  1  1-  i  »  I  I  L  !  l1 
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FIG.  Kji  ruuiii  icrnpcraiurc  saiuraiion 

torcsistancc  on  (a)  Co  interface  layer  thickness, 
sandwiches  of  the  form  Si/Py(53 -//)/Co(r/)/Cu(32), 

Py(22 -r/)/FeMn(90)/Cu(10),  (b)  distance  of  a  5  A  thic 
Co  layer  from  the  Py/Cu  interfaces  in  sandwiches  of  i\ 
form  Si/Py(49-rf)/Co(5)/Py(i/)/Cu(30)/Py(t/)/Co(5)/Py(] 
"t/)/FeMn(90)/Cu(10),  and  (c)  Py  interface  layer  thicknes 
ti,  in  sandwiches  of  the  form  Si/Co(57 -r/)/Py(/py)/Cu(24 
Py(//)/Co(29 ~ //)/FeMn(100)/Cu(l0).  Note  layer  thickness^ 
are  in  angstroms.  q, 


2.  Soatial  image  of  the  eigenstates  of  a  quantum  corral.  (Ai 


POWER  THRESHOLDS  AND  FREQUENCY  SHIFTS  IN  THE  EXCITATION 
OF  MAGNETOSTATIC  VOLUME  WAVE  SOLITONS  IN  DELAY  LINE 

CONFIGURATIONS 
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OUTLINE 


THE  CONDITIONS  FOR  THE  ONSET  OF  MODULATIONAL 
INSTABILITY  AND  EXCITATION  OF  SOLITONS  IN  MAGNETOSTATIC 
VOLUME  WAVE  DELAY  LINES  WILL  BE  PRESENTED. 

THE  FREQUENCY  AND  WAVEVECTOR  SHIFTS  INDUCED  BY  THE 
NONLINEAR  EFFECTS  AND  THE  POWER  THRESHOLDS  FOR  SOLITONS 
IN  REAL  DEVICE  STRUCTURES  WILL  BE  PREDICTED. 


MODULATIONAL  INSTABILITY 


INSTABILITY  IS  A  NONLINEAR  EFFECT  DUE  TO  BOTH, 
THREE-  AND  FOUR-MAGNONS  INTERACTIONS. 

MODULATIONAL  INSTABILITY  IS  DUE  TO  FOUR-MAGNON  SELF¬ 
INTERACTION  PROCESSES  ONLY. 

THE  ACTIVATION  OF  THE  MULTI-MAGNONS  DECAY  PROCESSES 
IN  CW  REGIME  IS  A  THRESHOLD  EFFECT  WHICH  RESULTS  IN  A 
MULTI-STABLE  RESPONSE.  THAT  BEHAVIOUR  CAN  BE  DETECTED  BY 
MEANS  OF  POWER  MEASUREMENTS  OR  FREQUENCY  SPECTRUM 
MEASUREMENTS,  THUS  EVIDENTIATING  THE  CREATION  OF 
SATELLITE  FREQUENCIES. 

THE  SATELLITES  LOCATION  DETERMINES  THE  KIND  OF 
MULTI-MAGNONS  INTERACTION,  AND  THEIR  PRESENCE  CRITICALLY 
DEPENDS  ON  THE  FREQUENCY  AND  ON  THE  POWER  LEVELS. 

A  FURTHER  PARAMETER  TO  BE  ACCOUNTED  IS  THE  LIGHTHILL 
CRITERION  SATISFACTION.  IN  FACT,  MODULATIONAL  INSTABILITY 
CAN  BE  DESCRIBED  BY  THE  NONLINEAR  SCHROEDINGER  EQUATION 
(NLSE)  AND  THE  RATIO  BETWEEN  THE  NONLINEAR  AND  THE 
DISPERSIVE  TERM  MUST  BE  NEGATIVE. 


GROUND  PLANES  CONTRIBUTION  TO  THE  DISPERSION 


THE  DISPERSION  RELATIONS  FOR  MSW  ARE  IMPLICIT,  AND 
FROM  THEIR  ANALYSIS  IT  TURNS  OUT  THAT  THE  STRONGEST 
CONTRIBUTION  TO  THE  GROUP  VELOCITY  TREND  IS  DUE  TO  THE 
GROUND  PLANES  CONFIGURATION. 

ONE  GROUND  PLANE  IS  NECESSARY  TO  GUARANTEE  THE 
ELECTRICAL  MATCHING  OF  THE  MICROSTRIP  LINE  USED  TO  LAUNCH 
THE  MICROWAVE  SIGNAL  INTO  THE  MAGNETIC  FILM,  WHILE  THE 
PRESENCE  OF  A  SECOND  GROUND  PLANE  MIGHT  BE  INTRODUCED  IN 
REAL,  CLOSED  DEVICE  STRUCTURES. 

THE  INTRODUCTION  OF  A  SECOND  GROUND  PLANE  CHANGES 
DRAMATICALLY  THE  DISPERSION  TREND,  AND  THE  EXSISTENCE  OF 
NONLINEAR  WAVES  FULLFILLING  THE  LIGHTHILL  CRITERION  IS 
ALLOWED  ONLY  IN  WELL  DEFINED  WAVEVECTOR  INTERVALS. 


ENVELOPE  SOLITONS 


ENVELOPE  SOLITONS  OF  NLSE  ARE  EXCITED  IN  A  LOW  LOSS 
MAGNETIC  FILM  WHEN  MICROWAVE  PULSES  ARE  LAUNCHED  IN  A 
DELAY  LINE  DEVICE,  WHERE  THE  MAGNETIC  FILM  IS  THE 
WAVEGUIDE  FOR  THE  PROPAGATION  OF  THE  ENVELOPE. 

THE  SATISFACTION  OF  THE  LIGHTHILL  CRITERION  IS  THE 
BASIC  REQUIREMENT  TO  IDENTIFY  THE  SOLITONIC  NATURE  OF  THE 
TRAVELING  ENVELOPES. 

IN  MAGNETOSTATIC  WAVE  (MSW)  DELAY  LINES,  THE 
NONLINEAR  TERM  MANTAINS  THE  SAME  ORDER  OF  MAGNITUDE  OVER 
THE  ENTIRE  EXCITATION  BAND  DEFINED  BY  THE  DC  MAGNETIC 
BIAS  FIELD  STRENGTH.  ON  THE  OTHER  HAND,  THE  DISPERSION 
RELATIONS  FOR  MSW'S  CRITICALLY  DEPEND  ON  THE  GEOMETRY  OF 
THE  PROBLEM  AND  NOT  ONLY  ON  THE  BIAS  FIELD  VALUE. 
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FREQUENCY  LIMITS  FOR  THE  EXCITATION  OF  MSW  SOLITONS  -  1 


A  MICROWAVE  PULSES  TRAIN  IS  REPRESENTED  IN  THE 
FREQUENCY  SPECTRUM  ACCOUNTING  FOR  THE  CARRIER  FREQUENCY 
fc,  THE  PERIOD  T  AND  THE  PULSEWIDTH  r  (OR  THE  DUTY  CYCLE 

D=r/T)  . 


THE  ENVELOPE  IN  THE  FREQUENCY  DOMAIN  IS  CENTERED  AT 
fc,  AND  THE  MAIN  LOBE  IS  LIMITED  BY  fc±l/r ,  WITH  THE 
FREQUENCY  COMPONENTS  WITHIN  THE  LOBE  SEPARATED  EACH  OTHER 
BY  1/T. 

TO  AVOID  THE  CUTTING  OF  SOME  FREQUENCY  COMPONENTS  OF 
THE  PULSE,  AT  LEAST  THE  MAIN  LOBE  MUST  LIE  ALL  WITHIN  THE 
MSW  BAND  EXCITATION.  IF  ft  AND  fn  ARE  THE  LOW  AND  HIGH 
FREQUENCY  LIMITS  FOR  THE  EXCITATION  OF  MSW'S,  IT  WILL  BE: 


ft  <  fc  -  1/t,  and  fH  >  fc  +  I/t- 


BY  IMPOSING  r  =  20  ns,  fc  MUST  BE  FAR  FROM  THE  FREQUENCY 
LIMITS  AT  LEAST  50  MHz. 


FREQUENCY  LIMITS  FOR  THE  EXCITATION  OF  MSW  SOLITONS  -  2 


GENERALLY  SPEAKING,  A  WIDEBAND  RESPONSE  IS  REQUIRED 
TO  CLEARLY  IDENTIFY  THE  PULSE  RESHAPING,  BECAUSE  A 
NARROWBAND  FILTER  INTRODUCES  ITSELF  A  SHAPE  MODULATION, 
THUS  MAKING  DIFFICULT  TO  RECOGNIZE  THE  SOLITONIC 
BEHAVIOUR  IN  THE  TIME  DOMAIN. 

THE  IDEAL  SITUATION  FOR  MSW  DEVICES  IS 
ASYMPTOTICALLY  REACHED  WHEN  WIDEBAND  MICROSTRIP 
TRANSDUCERS  ARE  COUPLED  TO  SUFFICIENTLY  THICK  MAGNETIC 
GARNET  FILMS  CHARACTERIZED  BY  A  SMALL  LINEWIDTH. 

BY  EXCITING  THE  FILM  AT  A  FREQUENCY  FAR  ENOUGH  FROM 
THE  MSW  BAND  LIMITS,  THE  PULSE  WILL  SUFFER  THE  DISPERSION 
CONTRIBUTION,  AND  THE  RESHAPING  OF  THE  PULSE  WILL  BE  DUE 
ONLY  TO  THE  SPREAD  INDUCED  BY  THE  INTRINSIC  DISPERSIVE 
BEHAVIOUR  OF  THE  MSW  DEVICE. 

AT  SUFFICIENTLY  HIGH  POWER  LEVELS,  THE  INTERPLAY 
BETWEEN  NONLINEARITY  AND  DISPERSION  WILL  CHANGE  AGAIN  THE 
PULSE  SHAPE  IF  A  SOLITON  IS  EXCITED.  A  PROPER  DESIGN  OF 
THE  DEVICE  EASILY  ALLOWES  FOR  THE  NECESSARY  WIDEBAND 


RESPONSE. 


ELECTRICAL  DESIGN  CRITERIA 


BY  INTRODUCING  THE  ELECTRICAL  DESIGN  OF  THE  DEVICE, 
THE  FILTER  BANDWIDTH  IS  NARROWER  THAN  THE  MSW  LIMITS, 
WHICH  ARE  ONLY  THE  PHYSICAL  LIMITS  FOR  THE  EXCITATION  OF 
MSW'S  DEFINED  BY  MEANS  OF  THE  BIAS  FIELD  STRENGTH. 

THE  REAL  BANDWIDTH  IS  OBTAINED  BY  USING  THE  CURRENT 
DENSITY  FLOWING  IN  THE  MICROSTRIP  TRANSDUCER  USED  TO 
EXCITE  THE  MAGNETIC  GARNET  FILM.  THE  CURRENT  DENSITY  J  IS 
A  FUNCTION  OF  THE  MICROSTRIP  WIDTH  Wm  AND  OF  THE  EXCITED 
WAVEVECTOR  k,  FOLLOWING  THE  RELATION: 


J  5=  sin(Wmk/2)/ (Wmk/2)  . 


THE  MAIN  LOBE  OF  THE  FILTERING  RESPONSE  IS  LIMITED  BY 
Wmk/2=7r,  THAT  IS  TO  A  MAXIMUM  ALLOWED  k-VALUE  OF  27r/Wm. 

COMMENT 

IT  REMAINS  THE  UNAVOIDABLE  LIMITATION  INTRODUCED  BY  THE 
NARROWBAND  FILTERING  RESPONSE  IN  THE  PULSE  RESHAPING, 
AND,  FROM  AN  EXPERIMENTAL  POINT  OF  VIEW,  THE  RISE  TIME 
AND  THE  FALL  TIME  OF  THE  REAL  PULSES  AT  THE  INPUT  OF  THE 


DEVICE . 


honliotiar  scheoedinger  equation  and  frequency  shifts  -I 

A  possible  approach  in  obtaining  a  NLSE  for  MSW'S  in 

PURE  DIPOLAR  REGIME  IS  THAT  OF  ZVEZDIN  AND  POPKOV  (1983) 
based  on  the  variation  OF  THE  DISPERSION  RELATION  WITH 
respect  to  THE  WAUEVECTOR,  THE  FREQUENCY  AND  THE  SQUARE 
modulus  of  the  MICROWAVE  MAGNETIZATION  AMPLITUDE 

normalized  by  the  SATURATION  VALUE,  IN  THE  APPROXIMATION 

OF  SMALL  SIGNALS- 

an  alternative  METHOD  IS  IN  CALCULATING 

VARIATION  DIRECTLY  ON  THE  FREQUENCY.  FROM  BOTH 
COMPUTATIONS,  THE  RESULT  IS  A  NLSE  WITH  THE  NONLINEAR  AND 
dispersive  terms  easily  obtainable  FROM  PARTIAL 

DERIVATIVES  OF  THE  DISPERSION  RELATION. 

.  A  PECULIAR  RESULT  IS  THAT  THE  TRAVELING  WAVE  IS 
affected  by  FREQUENCY  AND  WAVENUMBER  SHIFTS,  BOTH  RELATED 
TO  THE  NONLINEARITY  AND  THE  DISPERSION.  THAT  CORRESPONDS 
TO  A  MODULATION  OF  THE  CARRIER  FREQUENCY  AND  WAVENUMBER, 
WHOSE  VALUES  ARE  CHANGED  WITH  RESPECT  TO  THE  ORIGINAL 
ONES  IN  THE  LINEAR  REGIME. 


//X 


NONLINEAR  SCHROEDINGER  EQUATION  AND  FREQUENCY  SHIFTS  -2 


THE  CARRIER  WAVENUMBER  kc  WILL  BE  CHANGED  IN 


kc'  =  kc  +  K 


AND  THE  CARRIER  FREQUENCY  Wc  BECOMES 


Wc'  =  Wc  +  0, 


WHERE  0  IS  THE  FREQUENCY  SHIFT  AND  k  IS  THE  WAVENUMBER 
SHIFT,  AND  THE  CONDITIONS  k  «  kc  AND  Q  «  wc  ARE 


SATISFIED 


N0NLINE3UI  SCHROEDINGER  EQUATION  AND  FREQUENCY  SHIFTS  -3 


BY  FOLLOWING  THE  ANALYSIS  PERFORMED  BY  HASEGAWA  AND 
TAPPER!  IN  THE  CASE  OF  OPTICAL  FIBERS,  THE  FREQUENCY 
VARIATION  CAN  BE  WRITTEN  AS: 


WHERE  |v>l^(m/Ms)2  IS  THE  QUADRATIC  DEVIATION  FROM  THE 
EQUILIBRIUM  SATURATION  MAGNETIZATION  Ms-  in  IS  THE 
AMPLITUDE  OF  THE  MICROWAVE  MAGNETIZATION  IN  THE  PLANE 
NORMAL  TO  THE  EXTERNAL  DC  BIAS  FIELD  Ho/  AND  THE  QUADRATIC 
TERM  ARISES  FROM  THE  TAYLOR  EXPANSION  OF  THE 
MAGNETIZATION  COMPONENT  ALONG  THE  DIRECTION  PARALLEL  TO 
Ho*  LET'S  ASSUME  Z  THIS  COMPONENT,  IT  WILL  BE: 


NONLINEAR  SCHROEDINGER  EQUATION  AND  FREQUENCY  SHIFTS  -4 


SINCE  THE  THE  SECOND  ORDER  TERM  IS  THE  LOWEST  NOT 
NEGLIGIBLE  TO  BE  CONSIDERED  FOR  THE  NONLINEAR 
CONTRIBUTION,  THE  SAME  ORDER  HAS  TO  BE  USED  FOR  THE 
DISPERSION. 


WHERE  6a)Disp  AND  Swnl  ARE  THE  CONTRIBUTIONS  TO  THE  SHIFT  DUE 
TO  THE  DISPERSION  AND  THE  NONLINEARITY,  RESPECTIVELY. 
BOTH  ARE  EXPANDED  UP  TO  THE  SECOND  ORDER  AND  ARE  DEFINED 
BY; 


Us 


NONLINEAR  SCHROEDINGER  EQUATION  AND  FREQUENCY  SHIFTS  -5 

AN  EXSTIMATION  OF  THE  WAVENUMBER  MODULATION  IN  THE 
CASE  OF  MODULATIONAL  INSTABILITY  CAN  BE  GIVEN  BY  ASSUMING 
THAT  THE  INSTABILITY  GROWTH  RATE  /3  GROWS  UP  TO  A  MAXIMUM 
VALUE  DEFINED  BY: 


/c  =  Ko  =  y2|^|/|  (3Vak2)  I  =  y2lwr|/|  (aV3k2)  I 


WHERE  u)^7AH/2  IS  THE  RELAXATION  FREQUENCY,  7  IS  THE 
GYROMAGNETIC  RATIO  AND  AH  IS  THE  MAGNETIC  FULL  LINEWIDTH. 

THE  LOCATION  OF  THE  SATELLITES  WITH  RESPECT  TO  THE 
INPUT  FREQUENCY  Wc  CAN  BE  ROUGHLY  EXSTIMATED  BY  THE 
QUANTITY  Q  =  Vg/K  IN  A  FIRST  ORDER  APPROXIMATION,  WHERE  Vg 
IS  THE  GROUP  VELOCITY  AND  0  «  Wc- 
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WHEN  SOLITONS  ARE  CONSIDERED,  IT  IS  DIFFICULT  TO 
SEPARATE  THE  CONTRIBUTIONS  OF  THE  DISPERSION  AND  THE 
NONLINEARITY,  BECAUSE  THE  COMPONENTS  OF  THE  PERMEABILITY 
TENSOR  ARE  THEMSELVES  AFFECTED  BY  THE  NONLINEARITY,  AND 
THEY  MAKE  PART  IN  THE  DEFINITION  OF  THE  PARTIAL 
DERIVATIVES  USED  TO  CALCULATE  THE  DISPERSION  QUANTITIES. 


THE  GROUP  VELOCITY  IS  CHANGED  AND,  FOR  THAT  REASON, 
IT  IS  ONLY  A  ROUGH  APPROXIMATION  TO  ASSUME  THAT 
DISPERSION  AND  NONLINEARITY  PLAY  INDEPENDENT  ROLES  IN  THE 
DEFINITION  OF  0. 

-  LET'S  ASSUME  THE  VALIDITY  OF  THIS  APPROXIMATION,  WE 
WILL  CALCULATE  THE  iSwnl  VALUE  AS  A  FUNCTION  OF  THE  EXCITED 


WAVEVECTOR. 


DISPERSION,  NONLINEARITY  AND  LIGHTHILL  CRITERION 

THE  FREQUENCY  BEHAVIOUR  OF  ONE-  AND  TWO-GROUND 
PLANES  CONFIGURATIONS  FOR  VOLUME  MSW  DELAY  LINES  WILL  BE 
SIMULATED,  TO  ESTABLISH  THE  FREQUENCY  LIMITS  ALLOWED  FOR 
THE  SOLITONS  EXCITATION.  BOTH,  FORWARD  AND  BACKWARD 
VOLUME  WAVES  WILL  BE  CONSIDERED. 

FROM  THE  ANALYSIS  OF  THOSE  PICTURES,  IT  TURNS  OUT 
THAT  THE  POSSIBILITY  TO  FULLFILL  THE  LIGHTHILL  CRITERION 
FOR  THE  EXCITATION  OF  ENVELOPE  SOLITONS  IS  STRICTLY 
RELATED  TO  THE  GEOMETRY  OF  THE  DEVICE,  AND  MAINLY  TO  THE 
PRESENCE  OF  A  SECOND  GROUND  PLANE  AND  TO  THE  FILM 


THICKNESS  VALUE. 
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NONLINEAR  FREQUENCY  SHIFT 


—  THE  SHIFT  FOR  THE  ONE-GROUND  PLANE  FVW'S  IS  ALWAYS 
POSITIVE  FOR  ALL  THE  ALLOWED  k-VALUES ,  WHILE  POSSIBLE 
NEGATIVE  VALUES  DO  NOT  MAKE  SENSE,  BECAUSE  IN  THAT  CASE 
THE  LIGHTHILL  CRITERION  IS  NOT  FULLFILLED. 

TWO-GROUND  PLANE  FVW'S  ARE  CHARACTERIZED  BY  THE 
SATISFACTION  OF  THE  LIGHTHILL  CRITERION  ONLY  FOR  HIGH  k- 
VALUES  AND  THICK  FILMS.  FOR  THOSE  REASONS,  FVW  MSW'S 
SOLITONS  SEEM  LESS  FAVOURED  FOR  APPLICATIONS,  BECAUSE 
THEY  MEET  THE  RIGHT  SPECIFICATIONS  TOO  FAR  FROM  THE 
FREQUENCIES  WHERE  THE  DEVICE  USUALLY  EXHIBITS  THE  LOWEST 

INSERTION  LOSS. 

BVW  SOLITONS  ARE  LESS  SENSITIVE  TO  THE  FILM 
THICKNESS  AND  TO  THE  GROUND  PLANES .  THEY  HAVE  HAVE  ALWAYS 
NEGATIVE  SHIFTS,  UNLESS  IN  THOSE  REGIONS  WHERE  LIGHTHILL 
PROHIBITS  THEIR  CREATION .  ROUGHLY ,  THE  EXPECTED  SHIFT  FOR 
WAVEVECTORS  k  >  40  cm*^  IS  Q  =  -10^  sec^,  CORRESPONDING  TO 
FEW  HUNDREDS  OF  kHz. 

OWING  TO  THE  EXSISTENCE  OF  ALLOWED  AND  FORBIDDEN 
REGIONS,  THE  FREQUENCY  SPECTRUM  OF  A  PULSE  PROPAGATING  IN 
THE  DELAY  LINE  MUST  LIE  WITHIN  THE  ALLOWED  REGIONS. 
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THE  THRESHOLD  POWERS  FOR  THE  EXCITATION  OF  SOLITONS 
IN  A  MAGNETIC  GARNET  FILM  CAN  BE  EASILY  DERIVED  BY 
CALCULATING  THE  POYNTING  VECTOR  AND  BY  PERFORMING  AN 
INTEGRATION  OVER  THE  FILM  CROSS  SECTION. 

THIS  APPROACH  IS  RIGOROUSLY  VALID  BY  ASSUMING  AN 
HOMOGENOUS  MAGNETIZATION  OVER  THE  ENTIRE  FILM  WIDTH,  AND 
IT  LOSSES  MEANING  WHEN  NARROW  FILMS  ARE  CONSIDERED. 

IN  FACT,  THE  DEMAGNETIZATION  FACTOR  IS  ALMOST 
CONSTANT  ONLY  WHEN  THE  FILM  IS  WIDE  ENOUGH  WITH  RESPECT 
TO  ITS  THICKNESS.  FURTHERMORE,  A  NARROW  FILM  IS  AFFECTED 
BY  WIDTH  MODES,  WHICH  MAKE  WORST  THE  BAND  SHAPE  BY 
INTRODUCING  TRANSVERESE  RESONANCES  AND  THE  NOTCHES 
CORRESPONDING  TO  THEIR  FREQUENCY  VALUES. 

FOR  ABOVE  REASONS,  IN  MANY  PRACTICAL  APPLICATIONS  A 
FILM  WIDE  AT  LEAST  4  mm  MUST  BE  CONSIDERED  TO  ALLOW  FOR  A 
RELIABLE  PREDICTION  BY  USING  THE  POYNTING  APPROACH. 
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THRESHOLD  POWERS  -2 


THE  THRESHOLD  POWER  Pthrc  HAS  BEEN  COMPUTED  BY  USING 
THE  CLASSICAL  DEFINITION  FOUND  ALSO  IN  THE  ZVEZDIN  AND 
POPKOV  PAPER  (1983) . 

WE  HAVE  NOT  LIMITED  THE  CALCULATION  AT  k=0  (WHICH  IS 
NOT,  RIGOROUSLY  SPEAKING,  WITHIN  THE  MSW  APPROXIMATION), 
BUT  WE  HAVE  DERIVED  THE  k-DEPENDENCE  OF  Pthre  INSIDE  THE 
BAND  OF  THE  ALLOWED  k-VALUES .  THE  k=0  CONDITION  HAS  BEEN 
USED  ONLY  AS  A  CHECK  OF  THE  COMPARISON  BETWEEN  OUR 
RESULTS  AND  THOSE  OBTAINED  IN  THE  k=0  LIMIT  BY  ZVEZDIN 
AND  POPKOV. 


-  THE  CLASSICAL  COMPUTATION  IS  PERFORMED  BY  USING  THE 
FOLLOWING  FORMULA: 


Pthre  -  RS( 

p-dS  }= 

=  (o,/8n)  *  (-(1+Xl)/  (X1^X2^)  )  ■  (47rMs)2*d-W  \<p\^/^ 

WHERE  w  IS  THE  RADIAN  FREQUENCY,  xi  AND  X2  ARE  THE  DIAGONAL 
ELEMENTS  OF  THE  MAGNETIC  SUSCEPTIBILITY  TENSOR,  Ms  IS  THE 
SATURATION  MAGNETIZATION,  d  IS  THE  FILM  THICKNESS,  w  IS 
THE  FILM  WIDTH  AND  k  IS  THE  EXCITED  WAVEVECTOR. 


THRESHOLD  POWERS  -2A 


THE  THRESHOLD  AMPLITUDE  FOR  THE  NORMALIZED  MICROWAVE 
MAGNETIZATION  IS  GIVEN  BY: 


WHERE  Vg  IS  THE  GROUP  VELOCITY  AND  t  THE  PULSEWIDTH. 


IZS' 


IXi> 


wavevector  (cm-i) 


THRESHOLD  POWERS  -3 

FROM  THE  COMPUTATION  OF  THE  THRESHOLD  POWER  P,hrt  FOR 
THE  MSW  SOLITONS,  IT  TURNS  OUT  THAT: 

(i)  Pthre  DEPENDS  ON  THE  k-VALUE  AND  ON  THE  APPLIED  FIELD, 
AND  IT  IS  PARTICULARLY  EVIDENT  IN  THE  LIMIT  FOR  VANISHING 

k. 

(ii)  A  DRASTIC  CHANGE  IS  PREDICTED  IN  THE  TREND  AND  IN 
THE  ABSOLUTE  VALUE  OF  P,hre  PASSING  FROM  THE  NO  GROUND 
PLANES  CONDITION  TO  THE  ONE  GROUND  PLANE  GEOMETRY. 

(i)  IS  MAINLY  DUE  TO  THE  DISPERSION  TREND  AND  TO  THE 
ABSOLUTE  WORKING  FREQUENCY  DEPENDENCE  OF  Pthre-  THAT  IS 
ALSO  A  FACT  TO  BE  ACCOUNTED  FROM  AN  ENGINEERING  POINT  OF 
VIEW,  BECAUSE  THE  SAME  WAVEVECTOR  WILL  HAVE  NOT  THE  SAME 
THRESHOLD  AT  DIFFERENT  FREQUENCIES.  THIS  IS  ONLY  A 
VIRTUAL  PROBLEM  IN  THE  APPLICATION  OF  SOLITONS  TO  MSW 
DELAY  LINES,  BECAUSE  IT  IS  SUFFICIENT  TO  TAKE  INTO 
ACCOUNT  THE  Pthre  INCREASE  VS  THE  OPERATIVE  FREQUENCY. 

THE  (ii)  EFFECT  IS  TIGHTLY  RELATED  TO  THE  GEOMETRY 
OF  REAL  DEVICES,  AND  IT  HAS  TO  BE  CAREFULLY  ACCOUNTED  FOR 
APPLICATIONS . 
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TO  ACCOUNT  FOR  BOTH  THE  OUTLINED  CONTRIBUTIONS,  A 
DIFFERENT  Ho  VALUE  HAS  BEEN  IMPOSED  FOR  THE  SIMULATION, 
LEAVING  UNAFFECTED  THE  GEOMETRY.  THE  TREND  OF  Pthrc  vs  k  IS 
UNCHANGED,  BUT  A  FACTOR  CLOSE  TO  10  MAKES  THE  DIFFERENCE 
BETWEEN  THE  TWO  CURVES. 


).00  100.00  150.00  200.00 

wavevector  (cm-1) 


).00  100.00  150.00 

wavevector  (cm-l) 


THRESHOLD  POWERS  AND  ELECTRICAL  COUPLING 


THE  MICROSTRIP  TRANSDUCERS  USED  TO  EXCITE  THE  MSW 
DELAY  LINE  ARE 

1)  BIDIRECTIONAL 

2)  COUPLED  TO  MSW  UP  TO  A  MAXIMUM  ALLOWED 
EFFICIENCY. 

ENERGY  IS  LOST  IN  THE  COUPLING,  THUS  ENHANCING  THE 
INPUT  POWER  NECESSARY  TO  HAVE  AT  LEAST  THE  FIRST 
THRESHOLD  VALUE  IN  THE  FILM. 

FROM  THE  ELECTRICAL  QUANTITIES  INVOLVED  IN  THE 
ELECTRICAL  MATCHING,  A  COUPLING  COEFFICIENT  CAN  BE 
DERIVED  WHICH  MUST  BE  ACCOUNTED  IF  THE  REAL  THRESHOLD  FOR 
THE  DEVICE  HAS  TO  BE  ACCOUNTED. 

AS  AN  ORDER  OF  MAGNITUDE,  A  DELAY  LINE  LOSSING  -10 
dB  MUST  BE  EXCITED  BY  A  PULSE  TEN  TIMES  HIGHER  WITH 
RESPECT  TO  THE  THRESHOLD  POWER  NECESSARY  TO  CREATE  THE 


FIRST  SOLITON 


CONCLUSIONS 


THE  BASIC  CONDITIONS  FOR  DISPERSION  AND  NONLINEARITY 
TO  EXCITE  MAGNETOSTATIC  VOLUME  WAVE  SOLITONS  IN  DELAY 
LINE  CONFIGURATIONS  HAVE  BEEN  DISCUSSED.  PARTICULAR  CARE 
HAS  BEEN  DEDICATED  TO  THE  GROUND  PLANES  CONTRIBUTION. 

THE  MODULATION  IN  FREQUENCY  INDUCED  BY  THE  SOLITONS 
CREATION  HAS  BEEN  EVALUATED. 

THE  THRESHOLD  POWERS  FOR  BOTH,  FORWARD  AND  BACKWARD 
WAVES  HAVE  BEEN  DERIVED. 

THE  THRESHOLD  DEPENDS  ON  THE  GEOMETRY  AND  ON  THE 
ABSOLUTE  FREQUENCY. 

BVW'S  BEHAVIOUR  IS  NOT  CHANGED  BY  THE  PRESENCE  OF  A 
SECOND  GROUND  PLANE. 
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ABSTRACT 


We  report  the  observation  of  transient  chaotic  spin-wave  auto-oscillations  in  X- 
band  pulsed  subsidiary-resonance  experiments  in  a  YIG  sphere.  Chaotic  transients  to  a 
periodic  attractor  have  been  observed  near  a  critical  microwave  power  pc  ~  10  dB 
relative  to  the  Suhl  instability  threshold.  By  an  averaging  of  100  randomly  chosen 
initial  conditions  at  each  power,  we  estimate  a  critical  exponent  y  ~  0.57.  In  addition, 
we  demonstrate  that  the  results  can  be  qualitatively  understood  within  the  framework  of 
the  standard  two-mode  model.  Numerical  simulations  yield,  for  a  particular  set  of 
parameters,  chaotic  transients  obeying  the  predicted  scaling  law  with  a  critical  exponent 
y-  n.i4. 


1.  INTRODUCTION 


The  transition  to  sustained  chaotic  behavior  via  chaotic  transients  was 
discovered  by  Yorke  and  Yorke  in  the  Lorenz  model  in  the  late  70s.'  Several  years 
later,  Qrebogi,  Ott  and  Yorke^  identified  transient  chaos  as  a  particular  manifestation  of 
a  r,^t^]er  general  phenoinenon  that  occpr  in  dissipative  dynamical  systems,  the  so  called 
crk'is.  Three  types  of  crisis  have  been  characterized^,  namely,  attractor  destr-uction, 
attractor  met  gin  and  attractor  widening.  A  chaotic  attractor  is  destroyed  when,  as  a 
.systeni  parameter  p  reaches  a  critical  value p^. ,  it  collides  with  the  stable  manifold  of  an 
unstable  periodic  orbit.  The  characteristic  behavior  in  this  case  is  the  existence  of 
chaotic  transients  suddenly  followed  by  periodic  orbits.  The  duration  of  a  chaotic 
transient  depends  sensitively  on  the  initial  conditions  and  can  be  remarkably  long.  The 
average  lifetime  depends  upon  the  system  control  parameter  as  <  r  >  -  I  /?  -  Pc  I's 
where  y  is  the  critical  exponent  of  the  chaotic  transient. 

Finite-time  chaotic  behavior  has  been  observed  in  a  number  of  physical 
systems. 3  However,  only  a  few  experiments  have  been  analysed  in  detail.  As  far  as 
spin  waves  are  concerned,  the  only  results  available  on  such  transients  were  obtained  by 
Carroll,  Pecora  and  Rachford  **  a  few  years  ago,  in  subsidiary-resonance  experiments  in 
a  YIG  sphere.  At  driving  frequencies  between  2.0  and  3.4  GHz,  their  experimental 
findings  suggested  the  existence  of  "multiple  attractors",  while  numerical  simulations 
suggested  the  need  of  more  than  three  interacting  spin-wave  modes  to  explain  the  long- 
lived  transients  they  had  observed  in  the  experiments. 

In  this  paper  we  report  on  transient  chaos  in  YIG  spheres  at  higher  driving 
frequencies  (8.0  -  10.0  GHz)  at  room  temperature.  Pulsed  experiments  have  revealed 
the  existence  of  long-lived  chaotic  transients  involving  "single  attractors"  in  the  spin- 
wave  auto-oscillations  ("  1  MHz),  9  to  10  dB  beyond  the  Suhl  instability  threshold.  In 
addition,  we  present  numerical  results  of  a  two-spin-wave-mode  calculation  in  good 


qualitative  agreement  with  the  experiments.  By  an  averaging  of  N  randomly  chosen 
initial  conditions  at  each  power,  we  estimate  a  critical  exponent  in  the  range  0.5  <  y  < 
0.6  both  in  the  experiments  (N  =  100)  and  in  the  numerical  simulations  (N  =  400). 

2.  EXPERIMENTAL  RESULTS 

Our  experiments  were  carried  out  with  single-crystal  samples  of  yttrium  iron 
garnet  (YIG)  at  room  temperature.  Here  we  present  results  in  a  1-mm  YIG  sphere, 
which  is  placed  at  the  center  of  a  critically  coupled  rectangular  TE|Q2  microwave  cavity 
(Q  ~  2000).  An  electromagntet  provides  the  static  magnetic  field  Hq  perpendicular  to 
the  microwave  magnetic  field  h,  in  the  subsidiary  resonance  region.  The  power  is 
provided  by  a  10-W  traveling-wave-tube  (TWT)  amplifier  fed  by  a  solid-state  tunable 
YIG  oscillator.  The  frequency  {fp  =  8.9  GHz)  is  stabilized  by  an  external  crystal 
oscillator  and  manually  adjusted  to  the  center  of  the  cavity  resonance.  The  radiation 
power  is  controlled  with  a  variable  precision  attenuator  and  directed  by  a  circulator  to 
the  cavity,  where  it  drives  spin  waves  in  the  sample.  The  reflected  microwave  signal  is 
then  detected  with  a  sensitive  Schottkyrhanrier  diode  at  the  output  port  of  the  circulator 
and  recorded  at  intervals  of  0. 1  ps  using  a  commercial  digitizer.  In  order  to  study  the 
transient  response,  the  microwave  is  pulsed  by  a  p-i-n  modulator  placed  before  the 
TWT  amplifier.  Pulses  up  to  800  ps  long  were  comfortably  used  with  no  detectable 
heating  effects.  The  results  presented  here  were  observed  at  a  fixed  magnetic  field  Hq  = 
i8.”l0  Oe,  parallel  to  the  [100]  crystal  axis.  At  low-power  levels  the  pulse  reflected 
fitnn  the  cavity  hfis  essentially  the  same  shape  as  the  incoming  microwave  pulse.  As 
the  power  is  increased,  abrupt  changes  in  the  pulse  shape  occur  due  to  spin-wave 
instabilities  at  subsequent  thresholds,  namely,  the  Suhl  instability  threshold  (SIT),  , 
the  spin-wave  auto-oscillation  threshold  ,  ,  and  a  sequence  of  bifurcations  that  lead 

to  spin-wave  chaos.^  Henceforth  we  will  consider  R  =  hlh^  as  our  control  parameter. 


Well  above  the  SIT  {R  =  1),  we  have  observed  the  destruction  of  the  strange  attractor 
due  to  the  collision  between  the  chaotic  attractor  and  a  coexisting  unstable  periodic 
orbit.  This  collision  occurs  at  the  critical  value  =  2.88.  For  R>  R^^^  observe  an 
intermittent  signal  with  laminar  (periodic)  regions  interrupted  by  chaotic  bursts.  For 
2.82  <  R  <  Rc,  <he  steady-state  response  is  periodic  after  a  chaotic  transient  whose 
duration  swings  intermittently,  as  we  observe  on  the  screen  of  an  oscilloscope.  Figure  1 
shows  a  frozen  digitized  version  of  a  chaotic  transient  with  duration  T-  618  ps,  with  R 
=  2.84.  The  measured  average  duration  as  a  function  of  the  control  parameter  is  shown 
in  the  log-log  plot  of  Figure  2  (solid  circles).  From  a  linear  regression  (solid  line)  we 
estimate  a  critical  exponent  y  =  0.57. 


NUMERICAL  SIMULATIONS  AND  DISCUSSION 


The  efforts  in  modelling  the  instabilities  observed  in  high-power  ferromagnetic 
resonance  experiments  are  based  on  a  theory  introduced  by  Suhl^  in  the  mid  50s,  and 
have  been  phrased  jn  terms  of  the  excitation  and  the  nonlinear  interaction  of  spin  waves 
ill  the  sample.  The  number  and  the  nature  of  the  excited  modes  are  essentially  the 
hitherto  unsolved  problems  in  this  longstanding  subject.  We  have  recently  provided 
experimental  results  that  strongly  support  the  two-mode  model  (TMM),5  at  least  for  the 
interval  h  <  h  <  h~.  In  this  section  we  present  numerical  results  with  the  TMM  in 
good  qualitative  agreement  with  the  transient  experiments  described  above.  The 
particular  set  of  parameters  are  the  same  we  have  originally  used  to  explain  the  unusual 
spectra  observed  in  parallel-pumping  experiments.^  In  this  case,  the  model  gives  higher 
values  for  the  auto-oscillation  frequency  and  threshold,  but  shares  interesting  qualitative 
features  with  the  experiments.  We  refer  to  reference  [5]  for  a  detailed  description  on 
the  microscopic  approach  with  the  TMM.  Figure  3a  shows  a  chaotic  transient  in  the 
phase  plane  rt]  vs.  n2,  where  nj  and  ^2  are  the  populations  of  the  two  excited  spin-wave 
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modes,  for  R  =  8.10000506,  We  have  chosen  an  initial  condition  inside  the  strange 
attractor  for  which  the  solutions  are  attracted  during  a  time  interval  JiT-  87.7,  where 
Yl  is  the  phenomenological  relaxation  rate  of  mode  1.  The  steady-state  solution  is  the 
period-3  attractor  shown  in  Figure  3b  for  90.0  <  y,/  <  3000.0.  This  scenario  is  quite 
similar  to  the  experiments:  By  decreasing  R,  an  intermittent  (chaos  period.3) 
solution  is  destroyed  at  the  critical  value  R^  =  8.100055...  below  which  there  is  a  sharp 
period-3  window.  Within  this  window  the  chaotic  bursts  no  longer  exist  and  the  steady- 
state  solution  is  reached  after  a  chaotic  transient,  whose  duration  is  very  sensitive  to  the 
initial  conditions.  By  an  averaging  of  400  randomly  chosen  initial  conditions  within  a 
small  volufTte  on  the  strange  attractor  at  each  value  of  the  control  parameter,  we  have 
obtained  the  result  shown  in  Figure  4  (solid  circles).  The  straight  line  is  a  fit  from 
which  we  estimate  critical  exponent  y  =  0.54.  One  might  guess  at  first  that  we  are,  in 
both  experiment  and  simulation,  close  to  a  situation  where  the  system  could  be 
described  by  a  one-dimensional  map  with  a  quadratic  extremum,  since  in  this  case  y  has 
the  universal  value  0.5.  However,  we  stress  that  preliminary  studies  with  return  maps 

do  not  point  in  this  direction  and  further  investigations  are  under  way. 

This  work  has  been  supported  by  the  Brazilian  agencies  CNPq,  FINEP,  and 
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FIGURE  CAPTIONS 


FIG.  1:  Detected  microwave  absorption  vs.  time  in  a  YIG  sphere  at  subsidiary 
resonance  with  Hp=  1830  Oe,  R  =  2.84  andy^  =  8.9  GHz. 

F!G.  2-  Measured  average  duration  of  chaotic  transients  as  a  function  of  the  control 
parameter.  The  straight  line  is  a  fit  with  a  slope  -y  =  -0.57. 

FIG.  3:  Phase  plane  nj  vs.  n2  showing  numerical  solutions  of  the  two-mode  model,  as 
described  in  ref.  [7],  with  R  =  8.10000506:  (a)  transient  strange  attractor  for  0  <  Yit  < 
88.0  and  (b)  period-3  limit  cycle  for  90.0  <  yjt  <  3000.0. 

FIG  4:  Calculated  average  duration  of  chaotic  transients  as  a  function  of  the  control 
parameter.  The  straight  line  is  a  fit  with  a  slope  -y  =  -0.543. 
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abstract 

Microwave  pumped  spin-wave  instabilities  in  YIG  spheres  were  one  of  the  Brst 
experimental  systems  used  to  demonstrate  that  chaos  can  be  suppressed  by  small 
periodic  modulations  in  an  accessible  system  parameter.  Here  we  show  that  the 
equations  describing  two  interacting  spin-wave  modes  account  satisfactorily  for  the 
experimental  results,  as  long  as  the  field  modulation  is  appropriately  introduced  m 
the  model.  The  finite  detuning  parameters  expressing  the  boundary  conditions  of 
the  sample  provide  a  natural  way  for  introducing  the  field  modulation.  Thus  the 
present  results  constitute  additional  evidence  of  the  validity  of  the  two-mode  model 
with  momentum-non-conserving  driving  Hamiltonian  used  to  explain  the  sample  size 
dependence  of  the  self-oscillations. 


INTRODUCTION 


Microwave  pumped  spin-wave  instabilities  in  YIG  spheres  were  one  of  the  first 
experimental  systems  used  to  demonstrate^  that  chaos  can  be  suppressed  by  small 
time-dependent  variations  in  an  accessible  system  parameter.  This  method  of  control¬ 
ling  chaos  is  based  on  the  fact  that  a  chaotic  attractor  usually  has  embedded  within 
it  an  infinite  number  of  unstable  periodic  orbits.  As  shown  by  Ott,  Grebogi,  and 
Yorke^  (OGY)  it  is  possible  to  tailor  the  time-dependence  of  a  relevant  parameter  so 
as  to  stabilize  a  chosen  unstable  orbit.  In  systems  with  relatively  slow  dynamics,  the 
parameter  variation  necessary  to  stabilize  the  orbit  can  be  determined  at  each  cycle 
according  to  the  prescription  of  OGY  and  the  result  fedback  into  the  system^-A  In  the 
case  of  spin- wave  instabilities,  where  the  self-oscillation  frequencies  lie  in  the  lange  of 
lOOkHz  —  lOMHz,  a  closed-loop  feedback  experiment  is  difficult  to  implement^.  So 
the  stabilization  of  the  orbits  heis  been  achieved^  with  a  periodic  modulation  signal 
provided  by  an  independent  oscillator.  In  general  terms  the  parametric  pertuibation 
scheme  of  Lima  and  Pettini®  accounts  for  the  observed  suppression  of  chaos,  however 
there  is  no  specific  model  for  the  spin-wave  case. 

In  this  paper  we  show  that  the  control  of  chaos  in  the  spin-wave  experiments  can 
be  explained  satisfactory  by  the  two-mode  model  employed  to  describe  spin-wave 
nonlinear  dynamics.  The  field  modulation  used  in  experiments  to  control  the  chaotic 
state  is  naturally  introduced  in  the  model  through  the  detuning  parameters  express¬ 
ing  the  boundaxy  conditions  which  have  been  shown  to  account  for  the  sample  size 
dependence  of  the  self-oscillations’^’®. 


EXPERIMENTS 


The  experiments  have  been  previously  described^  m  detail  and  we  only  present 
here  some  of  its  main  features  and  results.  They  were  earned  out  with  a  polished 
sphere  (diameter  1.0mm)  of  the  prototype  ferromagnet  yttrium  iron  garnet  (YIG) 
at  room  temperature  pumped  with  X-band  radiation  in  the  “subsidiary-resonance” 
configuration.  The  experimental  arrangement  is  the  same  used  in  spin-wave  pumping 
experiments^’”  except  for  a  loop  placed  inside  the  microwave  cavity  to  allow  the  mod¬ 
ulation  of  the  sample  biasing  field  H  =  Ho  +  6H  cos{27rU)  over  a  broad  frequency 

range  0  -  lOMHz,  typically  with  SH/Hq  ~ 

The  usual  experiment  to  study  spin-wave  phenomena  is  done  with  fixed  values  of 
Ho  and  varying  microwave  power.  At  low  power  levels  the  steady-state  reflection 
from  the  critically  coupled  cavity  is  negligible.  As  the  microwave  driving  field  h  is 
increased,  an  abrupt  change  in  reflection  occurs  at  the  Suhl  threshold  ha,  due  to  the 
parametric  excitation  of  a  magnon  pair  with  frequency  Uk  ^  and  wave  vector 

k  and  —k.  The  value  of  k  is  determined  by  the  frequency  ijJp,  the  field  Hq,  and  the 
condition  for  minimum  threshold,  which  depends  on  the  pumping  configuration. 

The  suppresion  of  chaos  has  been  observed  at  many  points  of  the  h  x  Hq  diagram. 
The  results  presented  here  were  observed  after  the  system  has  been  driven  to  a  full 
chaotic  regime  with  h  =  1.1  Oe  and  Ho  =  1750  Oe,  with  the  field  applied  in  the  [110] 
crystal  direction.  Fig.l(a)  shows  the  power  spectrum  of  a  chaotic  auto-oscillation 
for  8H  =  0,  displaying  a  characteristic  broadband  feature.  By  increasing  SH,  the 
spectrum  becomes  progressively  cleaner,  with  sharp  lines  characteristic  of  a  periodic 
signal.  Fig. 1(b)  shows  the  result  for  6H  =  0.435  Oe  and  fm  —  1480  kHz.  The  spec 
trum  in  this  case  corresponds  to  a  quasiperiodic  signal  with  fundamental  frequencies 
_  740  kHz  and  fi  =  1975  kHz,  and  a  subharmonic  component  at  /o/2.  In  Fig.2 


we  show  the  variation  of  the  critical  amplitude  6H'  necessary  to  control  chaos  with 
the  modulation  frequency  fm-  Notice  that  8H'  has  minima  at  values  commensurate 
with  the  fundamental  frequency  /o,  i.e.,  fo/fm  =  P/?.  where  p  and  q  are  integers.  In 
addition,  for  the  fm  range  shown  in  Fig. 2,  the  minima  are  strikingly  ordered  according 
to  the  devil’s  staircase,  i.e.,  the  values  between  p/q  and  p'/q'  are  characterized  by  the 
ratio  {p  +  p')/{(I  +  ^')-  The  suppression  of  chaos  in  the  shaded  regions  of  Fig.2  was  also 
confirmed  by  the  behavior  of  the  information  dimension  Di  and  the  metric  entropy 
K  of  the  attractor.  As  h  increases  at  fixed  fm  a.nd  crosses  the  boundaries  in  the 
8H  X  fm  diagram,  Di  and  K  approach  the  values  1  and  0  respectively,  characteristic 

of  a  periodic  orbit^. 

THEORETICAL  MODEL  AND  RESULTS 

The  usual  two-mode  model®’^^  for  describing  the  nonlinear  dynamics  considers  that 
the  interacting  spin-wave  system  is  driven  by  a  uniform  microwave  field  in  an  infinite 
medium,  so  that  the  driving  Hamiltonian  n'{i)  preserves  momentum.  However,  in 
a  finite  medium,  H^t)  can  be  nonzero  for  pairs  k,-P  for  which  k  and  k'  differ 
by  an  amount  of  the  order  of  the  reciprocal  of  the  sample  size  L.  In  this  case,  it 
can  be  shown^  that  a  new  driving  term  is  introduced  in  the  spin-wave  equations, 
having  magnitude  proportional  to  a  factor  aAJt  which  depends  on  the  wave-vector 
mismatch  Ak  =  k  -k'.  For  instance,  for  two  neigbouring  standing  waves  along 
X  with  Afcr  =  vIL^,  a^k  attains  its  maximum  value  2/7r  ~  0.64.  In  an  infinite 
medium,  a^k  =  0.  With  this  new  momentum  non-conserving  driving  term,  the 
equations  of  motion  for  the  slowly  varying  spin-wave  variables  cj.  for  two  neighboring 


modes  become^ 

Cl  =  (71  +  iAuji)  Cl  -  thpi  (c*  + 

-i2{Sic\c\  +  Si2Cic\-\-2Ti2C2C2C^  ,  (1) 

C2  =  (72  +  C2  —  ihp2  (c2  +  ae“'^/^ci) 

—  22  ^S'2C2C2  +  5i2C2C|  +  2Ti2CiCiC2j  ,  (2) 

where  Awi..  =  a;*,.  —  u;p/2,  a  =  ttAt  and  (i  =  /?i2  is  the  ph<ise  difference  between  modes 
1  and  2,  7*;  is  the  relaxation  rate,  pk  the  coupling  factor  between  mode  k  and  the 
pumping  field  h  and  the  5,  and  T,  are  the  nonlinear  coupling  parameters^. 

The  presence  of  one  mode  in  the  pumping  term  of  the  other  represents  a  modulation 
with  the  frequency  difference  Aw  =  Wi  —  Wj,  leading  to  self-oscillations  with  frequency 
that  depends  not  only  on  the  relaxation  rate  but  also  on  Aw.  For  two  neighboring 
modes  in  k  space,  Aw  oi  rfL,  which  accounts  for  the  observed  size  dependence  of 
the  spin-wave  self-oscillations.  When  a  small  ac  field  variation  is  superimposed  to  the 
static  field,  the  spin-wave  detuning  parameters  become 

Awfc  =  wfc  -  Y  -h  'y6H  cos{2Trfmt)  .  (3) 

Equations  (1)  and  (2)  with  the  time-dependent  detunings  (3)  constitute  our  model 
to  explain  the  suppression  of  chaos.  Those  equations  have  been  integrated  numerically 
in  a  SUN  Sparc  2  workstation  with  a  Runge-Kutta  subroutine.  Chaotic  dynamics  and 
control  of  chaos  with  field  modulation  are  observed  with  many  sets  of  parameters. 
The  results  described  below  were  obtained  with  71  =  72,  Pi  =  P2,  Awi/71  =  0.2, 
AW2/71  =  -0.5,  5i/7i  =  82/^1  =  0.5,  Su/'fi  =  -0.2,  Tn/^i  =  -0.4,  a  =  0.65 
and  j3  =  TT.  With  these  parameters,  auto-oscillations  with  frequency  /o  —  0.497i 
(with  no  field  modulation,  8H  =  0)  develop  at  a  normalized  microwave  amplitude 
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^  ^  =  1.32  (the  Suhl  threshold  \s  R  =  1.0)  and  chaos  sets  in  at  i?  -  1.S08. 

Fig. 3(a)  shows  the  power  spectrum  of  mode  1  amplitude,  rii  =  CiCi,  for  R  —  1.83 
and  5H  =  0,  characterizing  a  chaotic  state.  With  the  modulation  turned  on  at  some 
frequency  values,  if  SH  exceeds  a  critical  amplitude  6H\  chaos  is  suppressed  after  a 
transient  time.  Fig.3(b)  shows  the  power  spectrum  for  =  O.SOqi  and  SH  =  O.471/7 
(7  is  the  gyromagnetic  ratio),  demonstrating  that  the  signal  has  become  periodic 
with  frequency  0.56  71-  The  boundaries  of  the  stability  regions  in  the  modulation- 
amplitude  X  frequency  plane,  shown  in  Fig.4,  display  tongues  like  the  experimental 

results  of  Fig. 2. 

The  results  shown  in  Figs.  3  and  4  were  obtained  by  integrating  the  equations  of 
motion  with  the  field  modulation  on  for  a  length  of  time  on  the  order  of  200  cycles 
of  the  oscillation,  so  that  the  response  has  reached  steady-state.  .Actually,  the  stable 
orbit  resulting  from  the  field  modulation  is  preceded  by  a  chaotic  transient  in  which 
orbits  similar  to  those  of  the  uncontrolled  chaotic  attractor  dominate^.  The  duration 
T  of  such  chaotic  transient  depends  sensitively  on  the  initial  conditions.  For  a  given 
modulation,  we  calculate  the  average  duration  <  r  >  for  256  randomly  chosen  initial 
conditions.  For  =  0.56  71  this  average  time  is  found  to  scale  with  the  modulation 
amplitude  as  <  r  >~  {SH-SHT'^.  where  7  ^  0.287.  A  scaling  relation  of  this  type 
was  predicted  by  OGY,  though  with  a  different  value  of  the  critical  exponent. 


CONCLUSIONS 

The  experimental  observations  on  controlling  chaos  in  spin-wave  instabilities  driven 
by  microwave  radiation  in  YIG  spheres,  by  means  of  a  small  periodic  variation  in  the 
applied  magnetic  field,  are  explained  by  a  two-mode  model.  The  field  modulation 
is  introduced  in  the  model  through  the  same  detuning  parameters  used  to  explain 
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the  sample  size  dependence  of  the  self-oscillations.  Hence,  the  present  results  pro 
vide  additional  evidence  of  the  validity  of  the  two-mode  model  with  momentum-non 
conserving  driving  Hamiltonian  for  describing  spin-wave  dynamics. 
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FIGURE  CAPTIONS 


FIG.  1.  Power  spectra  of  observed  auto-oscillations  for  different  values  of  the  amplitude 
6H  of  the  field  modulation,  (a)  Chaos  for  SH  =  0.  (b)  Chaos  under  control  with  6H  = 
0.435  Oe  and  /m  =  1480  kHz  [Ref.l]. 


FIG.  2.  Critical  modulation  amplitude  SH*  vs.  modulation  frequency  fm-  The  bound 
aries  between  the  chaotic  and  controlled  regions  have  minima  at  values  of  commensurate 
with  the  fundamental  frequency  /o  at  ratios  indicated  at  the  top  [Ref.l]. 

FIG.  3.  Fourier  transform  of  mode  1  amplitude  without  (a)  and  with  field  modulation 
(b),  SH  =  0.4  71/7,  fm  =  0.56  71.  The  clean  spectrum  in  (b)  demonstrates  suppression  of 

chaos  with  the  two-mode  model. 


FIG.  4.  Critical  amplitude  SH*  in  units  of  71/7  obtained  with  the  model  and  parameters 


described  in  the  text. 


SAMPLE-SIZE  EFFECT 
IN  SPIN  WAVE 
AUTO-OSCILLATIONS 


A.N.  Slavin 

Department  of  Physics,  Oakland  University, 
Rochester,  MI  48309,  USA 


ABSTRACT 


The  theoretical  model  of  spatially  inhomogeneous 
unstable  collective  oscillations  in  a  system  of  spin  waves 
parametrically  excited  in  a  finite  magnetic  sample  is  used 
to  explain  the  sample-size  effect  in  spin  wave  auto¬ 
oscillations  under  parallel  pumping  observed  in 
yttrium-iron  garnet  (YIG)  spheres  and  films. 

It  is  shown  that  in  the  axially  symmetric  case  (YIG 
sphere)  the  model  demonstrates  good  quantitative 
agreement  with  the  experiment  for  both  auto-oscillation 
threshold  and  frequency. 

In  the  absence  of  axial  symmetry  (tangentially 
magnetized  YIG  film)  the  instability  criterion  for 
collective  oscillations  is  fulfilled  only  for  sufficiently 
small  values  of  the  bias  magnetic  field  H  (large  values  of 
parametric  spin  wave  wave  number  k  ).  In  this  case  the 
model  gives  only  qualitative  description  of  the  size 
dependence  of  the  auto-oscillation  threshold,  but 
correctly  describes  the  size  dependence  of  the  auto¬ 
oscillation  frequency. 


1.  Introduction 

Spin  wave  auto-oscillations  in  ferrites  above  the 
threshold  of  parametric  excitation  discovered  in  [1]  were 
studied  intensively  during  the  last  decade  mostly  because 
of  interesting  chaotic  dynamics  they  demonstrate  in  a 
strongly  nonlinear  regime  (see  e.g.  [2,  3,  4,  5]. 

The  experiments  performed  by  Rezende  et  al.  [6,  7] 
demonstrate  a  pronounced  sample-size  dependence  of 
auto-oscillation  frequency  and  threshold  that  was  not 
explained  by  the  traditional  theory  of  spin  wave  auto¬ 
oscillations  [8,  9]. 

In  our  previous  work  [10,  11]  we  have  developed  a 
model  of  spin  wave  auto-oscillations  under  parallel 
pumping  in  finite-size  ferromagnetic  samples  taking  into 
account  the  boundary  conditions  for  the  envelope  of  the 
packet  of  parametrically  excited  spin  waves  at  the 
boundaries  of  a  magnetic  sample.  This  model  explains  the 
sample-size  dependence  of  auto-oscillation  frequencies 
and  thresholds  observed  YIG  spheres  in  [6,  7]. 

In  our  present  paper  we  make  a  detailed 
comparison  of  this  model  with  the  experiments 
performed  by  Rezende  et  al.  in  YIG  spheres  [6,  7], 

describe  our  own  experiments  on  sample-size 
dependence  of  auto-oscillations  performed  in 

tangentially  magnetized  rectangular  YIG  film  samples, 
and  discuss  the  applications  of  our  theoretical  model  to 

the  case  of  film  samples. 
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2.  Theoretical  model 


The  theoretical  model  [10,  11]  that  we  are  using  to 
describe  the  sample-size  effect  in  spin  wave  auto¬ 
oscillations  is  an  extension  of  the  S-theory  of  auto¬ 
oscillations  in  unlimited  ferromagnetic  samples  [8].  In  the 
framework  of  this  theory  spin  wave  auto-oscillations 
are  interpreted  as  unstable  collective  oscillations  in  the 
system  of  parametrically  excited  spin  waves.  These 
oscillations  are  originated  from  the  four-wave  interaction 
between  parametric  spin  waves.  The  criterion  of  instability 
of  collective  oscillations  (i.e.  appearance  of  auto¬ 
oscillations)  can  be  formulated  as  (see  [8]) 

2R+1<0,  (1) 

where  R  =  T/S  ,  and  S  and  T  are  the  coefficients  of 
four-wave  interaction  between  spin  waves  . 

In  the  case  of  a  finite  magnetic  sample  we  assume  that 
the  collective  oscillations  in  the  system  of  parametrically 
excited  spin  waves  can  be  spatially  inhomogeneous  (i.e. 
can  propagate  in  the  sample)  and  the  amplitudes  Ak  of 
these  oscillations  (or  secondary  waves)  must  satisfy  certain 
boundary  conditions  at  the  sample  boundaries. 

For  simplicity  we  assume  that  only  one  size  of  the 
sample  is  relevant  -  the  size  D  in  the  direction  of 
propagation  of  parametric  spin  waves  excited  at  the 
threshold. 


D 


In  the  case  when  the  mean  free  path  1  -  vg/yk  of 
spinwaves  in  a  sample  is  much  smaller  than  the  sample 
size  D  (  1  «  D)  the  boundary  conditions  for  the 
oscillation  amplitudes  must  be  the  conditions  of 
"no  reflection"  at  the  boundary  (  Ak  =  0  ).  In  this 

particular  case  the  model  [10,11]  yields  the  following 
expressions  for  the  auto-oscillation  threshold  hose 
the  auto-oscillation  frequency  ®osc  threshold 

2  \Dyt 

Yk  \2Ji  +  l 

where  hth  is  the  threshold  of  parametric  excitation  of  • 

spin  waves  ,  Yk  is  the  relaxation  parameter,  and 
Vg  is  the  group  velocity  of  parametric  spin  waves  . 


The  four-wave  interaction  coefficients  T  and  S  (  and 
their  ratio)  can  be  calculated  from  the  general  expressions 
presented  in  [12].  We  note,  that  both  values  and  signs 
of  the  four-wave  interaction  coefficients  S  and  T 
depend  strongly  on  the  crystallographic  orientation  of 
the  sample  even  in  cubic  crystals  like  YIG,  where  spin 
wave  spectrum  is  practically  isotropic.  This  property  of  T 
and  S  the  giant  crystallographic  anisotropy  of  auto¬ 
oscillations  observed  in  experiments. 

Another.. interesting  feature  of  the  expressions  (2), 
(3)  is  that  in  the  axially  symmetric  case  (sphere)  the 
ratio  of  these  two  expressions  is  independent  of  both 
the  sample  size  D  and  the  bias  magnetic  field  H. 
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The  model  Eqs.{2),  (3)  is  applicable  when  parametric 
spin  waves  with  wave  numbers  k  >  0  ( H  <  He)  are 
excited  in  the  sample  and  describes  all  the  characteristic 
features  of  the  spin  wave  auto-oscillations  under  parallel 
pumping  (including  the  sample-size  effect). 


3.  Sample-Size  Effect  In  Spheres 

The  most  complete  experimental  investigation  of 
the  sample-size  effect  in  spin  wave  auto-oscillations 
under  parallel  pumping  in  YIG  spheres  has  been  reported 
by  Rezende  et  al.  in  [6].  The  results  for  the  threshold  of 
parametric  excitation  of  spin  waves  hth  ,  the  threshold  of 
spin  wave  auto-oscillations  hose  ?  the  auto¬ 
oscillation  frequency  fosc  ( fosc  =  COosc/^TT)  as  functions 
of  a  bias  magnetic  field  H  are  presented  in  Fig.  7  of 
Ref.  [6]  for  YIG  spheres  of  diameters  D  =  1  mm  and  0.52 
mm. 


Using  expressions  for  T  and  S  from  [12]  we  calculated 
the  ratio  G  Eq.(4)  for  the  conditions  of  the  experiments 
[6].  The  parameter  R  =  T/S  in  this  case  was  equal  to  R 
=  -  1.42,  while  the  ratio  G  turned  out  to  be  equal  to 


G  =  0.3  .  We  also  calculated  the  same  ratio  G  from  the 
experimental  data  presented  in  [6].  .  The  theoretical  and 
experimental  values  of  G  are  presented  in  Fig.  1. 

Tt  is  clear  from  Fig.l  that  although  the  numerical 
values  of  theoretical  and  experimental  ratio  G  are  differ¬ 
ent,  the  behavior  of  G  with  variation  of  D  and  H  is  simi¬ 
lar  :  in  both  theory  and  experiment  G  does  not  depend 
on  either  the  sample  size  D  or  the  bias  field  H. 

Using  Fig.l  we  determined  the  value  of  the  parameter 
R  for  which  the  theoretical  value  of  the  ratio  G  is  equal 
to  the  mean  average  of  experimental  values  G  =  0.8  . 
This  value  turned  out  to  be  equal  to  R  =  -  0.765.  We 
note,  that  the  instability  criterion  (1)  is  still  fulfilled  for 

this  value  of  R. 

Below  we  used  R  =  -  0.765  in  Eqs.(2),  (3)  to  calcu¬ 
late  the  theoretical  curves  hosc(H  )  ODosc(H)  •  The 
results  of  comparison  of  these  calculations  with  the 
experimental  data  from  [6]  are  presented  in  Fig.2. 

It  is  clear  from  Fig.2  that  our  model  Eqs.(2),  (3) 
demonstrates  good  quantitative  agreement  with  the  results 
of  the  experiment  [6]  and  explains  the  sample-size  effect 
in  spin  wave  auto-  oscillations  under  parallel  pumping  in 
YIG  spheres. 
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Fig.2  Comparison  of  the  results  of  experimental 
measurements  of  auto-oscillation  threshold  and  frequency 
performed  on  two  YIG  spheres  (D  =  1  mm  and  0.52  mm) 
in  [6]  with  the  results  of  theoretical  calculations  of  these 
characteristics  done  using  the  model  Eqs.(2),  (3)  for 
R  =  T/S  =  -  0.765. 


3.  Sample-Size  Effect  in  Films 

To  make  a  further  check  of  our  theoretical  model 
we  performed  a  series  of  experiments  on  rectangular  YIG 
film  samples  magnetized  in  the  film  plane  along  the  sample 
width  (W  =  2  mm).  The  idea  of  the  experiment  was  to 
change  the  length  of  the  sample  (D  =  0.8  -  6  mm)  and  to 
see  how  it  affects  the  frequency  of  spin  wave  auto¬ 
oscillations  COosc-  The  samples  were  cut  from  an 
epitaxial  YIG  film  (thickness  L  =0.11  mm)  of  the 
orientation  <1 1 1>  to  get  the  minimum  crystallographic 
anisotropy  in  the  film  plane.  The  spin  waves  were  excited 
by  parallel  pumping  in  a  TE  reflection  type  cavity  at  the 

frequency  CO  72711=  9.24  GHz  [14,  15].  The  microwave 
pumping  field  was  parallel  to  the  bias  magnetic  field  (and 
to  the  width  W  of  the  sample). 


The  experimental  results  for  tangentially  magnetized 
YIG  film  samples  of  different  lengths  D  are  presented  in 
Fig.3  a-d.  It  is  easy  to  see,  that,  contrary  to  the  case  of 
spherical  samples,  in  film  samples  the  sample-size  effect 
is  observed  only  for  sufficiently  small  values  of  the  bias 
magnetic  field  H  <  700  Oe  (or  for  sufficiently  large 
values  of  spin  wave  wavenumber  k>3  10^  cm’^  ). 
Theoretical  interpretation  of  the  size  effect  in  films  is 
more  complicated  than  in  spheres  for  the  following  reasons 

(i)  in  films  we  have  two  effective  sample  sizes  (D  and  L); 

(ii)  distribution  of  parametric  spin  waves  in  the  film  is 
not  axially  symmetric,  so  the  coefficients  S  and  T  are 
dependent  on  the  bias  field  H  . 

To  apply  our  one-dimensional  model  to  the  case  of  a 
tangentially  magnetized  film  sample  we  made  a  naive 
assumption  that  the  auto-oscillation  frequency  in  this 
sample  can  be  represented  in  the  form 

where  the  first  term  depends  on  the  sample  length  D  and 
is  given  by  Eq.(3),  while  the  second  term  is  determined 
by  the  film  thickness  L  and  by  small  effects  like  scattering 
of  auto-oscillations  on  each  other  and  on  random 
inhomogeneities  of  the  medium.  These  effects  are  described 
in  detail  in  [16] . 


D  = 


3  ^ 


Fig.3  Experimental  characteristics  of  spin  wave  auto¬ 
oscillations  in  tangentially  magnetized  YIG  film  samples  , 
(thickness  L  =0.11  mm,  width  W  =  2  mm,  orientation 
(111))  of  different  lengths  0.8  mm  <  D  <  6  mm  ; 

(a)  threshold  of  parametric  excitation  ;  (b)  threshold  of  ^ 

spin  wave  auto-oscillations  ;  (c)  firequency  of  spin  wave 
auto-oscillations;  (d)  frequency  of  spin  wave  auto- 
oscillations  for  several  intermediate  sample  lengths  D  m 
the  bias  field  region  (H  <  700  Oe)  where  the  sampk-^ze  • 


H  (Oe) 
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Fig.4  Comparison  of  theoretical  and  experimental  results 
for  frequency  of  spin  wave  auto-oscillations  in  tangentially 
magnetized  YIG  film  samples: 

(a)  dependence  of  the  frequency  on  the  bias  field  H  , 

(b)  dependence  of  the  frequency  on  the  sample  length  D 

/?5^ 


We  did  not  calculate  Q(L)  here  and  assumed  that 
this  term  is  equal  to  the  auto-oscillation  frequency 
190  kHz  in  the  sample  having  the  largest  length  D  =  6 
mm  (see  Fig.3  c,d). 

It  is  clear  from  Fig.3  b,  c  that  the  sample-size 
effect  in  film  samples  manifest  itself  only  for  H  <  700  Oe  . 
(The  nature  of  auto-oscillations  observed  in  the  bias  field 
region  H  >  700  Oe  will  be  discussed  elsewhere  .) 

To  explain  this  fact  we  assumed  that  for  H  >  700  Oe 
the  instability  criterion  Eq.(l)  is  broken  and  the 
coefficient  S(H  )  changes  its  sign  from  positive  to 
negative  at  H  =  700  Oe.  Using  this  assumption  we 
calculated  theoretical  dependencies  of  auto-oscillation 

frequency  fosc  =  O3osc/27r  on  the  sample  length  D  (at  a 
fixed  H  =  300  Oe)  and  on  the  bias  field  H  (at  two  fixed 
lengths  D  =  0.8  mm  and  D  =  2  mm)  from  Eqs.(3),  (5). 

The  results  of  comparison  of  these  calculations  with 
the  experimental  results  from  Fig.3  c,  d  are  presented  in 
Fig.4  a,b  .  It  is  clear  that  Eqs.  (3),  (5)  give  good 
qualitative  explanation  of  the  sample-size  effect  on  auto¬ 
oscillations  in  films. 


CONCLUSION 


The  above  presented  comparison  of  theoretical  and 
experimental  results  for  the  characteristics  of  spin 
wave  auto-oscillations  demonstrate  that  our  simple 
one-dimensional  theoretical  model  of  spin  wave  auto¬ 
oscillations  in  finite-size  magnetic  samples  works 
well  in  spheres,  but  gives  only  limited  qualitative 
description  of  the  sample-size  effect  on  spin  wave 
auto-oscillations  observed  in  tangentially  magnetized 
YIG  films.  Presently,  efforts  are  in  progress  to 
develop  a  two-dimensional  theoretical  model  of  spin 
wave  auto-oscillations  in  magnetic  films  capable  to 
give  the  quantitative  description  of  the  above 
reported  sample-size  effects. 
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Our  experiments  at  VHF/UHF  frequencies  and  at  low 
applied  fields  with  Soohoo  Associates  Magnetic  Domain 
Spectrometer  103  showed  that  Kittel's  resonance  equation 
should  be  modified  with  the  addition  of  to  read: 

to  =)Y|j  [  H^+H^+(N^-N^)Mg]  [  H^+H.  +  (Ny-Nz)Mg]  (D 

where  the  applied  field  H  may  include  anisotropy  fields. 
In  Eq.  (1),  N^,  Ny  and  are  demagnetizing  factors  in  the 
x,y,z  directions  and  M  is  the  saturation  magnetization. 

We  found  that  is  a  fraction  of  an  oersteds  for 
permalloy  films  and  10 's  of  oersteds  for  ferrite  discs. 

The  origin  of  has  been  attributed  by  E.  Schlomann 
to  the  directional  distribution  of  anisotropy  axes  in  a 
polycrystal  and  to  sample  porosity.^  Our  experimental 
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Modified  Kittel's  Resonance  Formula. 

Experimental  Setup  -  Soohoo  Associates'  Magnetic  Domain 
Spectrometer  System  MDSS103. 

Sample  Absorption  vs.  Applied  Field  Parallel  and 
Perpendicular  to  Easy  Axis . 

M-H  Loops. 

Absorption  vs.  Easy-Axis  Field  for  Ni/Fe  Film. 
Single-Pass  Barkhausen  Transitions  for  Ni/Fe  Film. 
Multiple-Pass  Barkhausen  Transitions  for  Ni/Fe  Film. 
Absorption  vs.  Hard-Axis  Field  for  Ni/Fe  Film. 

Easy-Axis  M— H  Loop  for  Ni/Fe  Film. 

Hard-Axis  M-H  Loop  for  Ni/Fe  Film. 

Hard-Axis  M-H  Loop  for  Ni/Fe  Film 
(without  domain-wall  motion). 

Internal  Field  for  Ni/Fe  Film  t2. 

Internal  Field  for  Ni/Fe  Film  Si. 

Absorption  vs .  Applied  Field  for  Polycrystalline 
Ferrite  Disc. 

Barkhausen  Transitions  for  Ferrite  Disc. 

Internal  Field  H^  for  Ferrite  Disc. 

Origin  of  Internal  Field. 

Resonance  Field  for  Different  Crystalgraphic 
Orientations . 

Internal  Fields  Expressions  due  to  Random  Crystallite 
Orientation  and  Porosity. 
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results  show  that  there  are  additional  contributions  to 

when  the  sample  is' unsaturated .  Whereas  our  experiments  at 

VHF/UHF  frequencies  are  performed  at  low  fields  of  0  to 

100  oe  which  encompass  non-saturation  and  changes  in 

2  3 

magnetic  states,  other  experiments  '  were  typically- 
performed  at  x-band  ('vlOGHz)  and  above  where  resonance 
fields  were  in  the  1000 's  of  oersteds  corresponding  to 
magnetic  saturation.  For  of  a  fraction  of  an  oersted 
that  we  measured,  it  would  not  be  possible  to  discern  at 
these  high  fields. 
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CRITICAL  BEHAVIOR  OF  MAGNETIC 

DOMAINS 

H.  SUHL 


University  of  California,  San  Diego 


Ordered  systems,  with  a  continuous  degree  of  symmetry,  invariably  exhibit  a 
singular  response  to  a  perturbation  coupling  to  the  degree  of  freedom  with 
that  symmetry.  (Goldstone  mode.)  A  very  simple  case  is  a  ferromagnet 
without  any  kind  of  anisotropy:  The  direction  of  its  magnetization  can  be 
changed  continuously  without  changing  the  energy  of  the  system. 
Correspondingly  the  transverse  susceptibility  at  low  field  diverges  as  the 
field  tends  to  zero.  In  symbols: 


H=H^+h,  M=M  +m 
To  lowest  order  in  h 


m  -h  I 


n 


0 


,  dmldh=ll 


H. 


0 


which  goes  to  infinity  as  the  biasing  field  goes  to  zero.  A  slightly  more 
elaborate  argument  shows  that  the  longitudinal  susceptibility  likewise 

diverges. 


II.  In  the  case  of  a  magnetic  sample  with  domain  walls  (i.e.  a  sample  in  a 
subsaturating  external  field)  there  is  likewise  “very  nearly”  a  continuous 
symmetry:  the  positions  of  the  domain  walls.  Neglecting  interaction 
between  them,  as  well  as  pinning,  they  can  be  moved  without  changing  the 
energy.  The  magnetization  component  along  the  field,  on  the  other  hand, 
can  change  as  the  result  of  wall  motion.  The  response  is  usually  quoted  as  a 
response  to  an  internal  magnetic  field.  In  an  obvious  notation: 


I  ^  ^ 

Xe  =  dM! dH^ = {dM!  dHMH  /  dH^) 

Xi=Xe>0--^zXe^ 

So  if  the  measured  Xq  equals  1/ then  X-  is  infinite 


I 


HI.  The  following  is  a  sketch  of  a  possible  theory  that  should  lead  to  the  correct 
critical  index  a  in  the  observed  law  « i  /  /»“. 


Hamiltonian: 

mm 

Partition  function : 

Z= 

Constrained  Partition  function  well  below  : 

where  Fi^=-Hg-M+^MNM-k^TSiA^ 

N  and  S  depend  only  weakly  on  T  and 
and  an  expansion  of  S  starts  with  terms  of 

fourth  order  in  the  components  of  M. 

The  first  two  terms  give  the  infinite  susceptibility  below  saturation,  but  the 
last  term  should  give  a  finite  result  « i  /  at  small  but  finite  h. 


2,2.4- 


Ilia.  The  unrenormalized  result  might  be  as  follows. 


Cylindrical  Symmetry  around  an  external  field 
in  the  z-direction. 

Mean  field  theory:  dF I  dM^=0,  that  is 


He =N^Mz-  kgTdS  /  dM^ 
so  if  S  is  neglected, 

If  leading  term  in  S  is 


H.=He-NzMz,  we  get 

I 


Mz=iH.I{k^Ta)^^,  X. 


2/3 


V. 


Initial 

state 


self  organized  crilicaj  state 


ground  elate 


SANDFILE 


-::iAD-ON  DOMAIN  WALL 


Formation  of  zig-zag  domain  wall  in  a  film  of  particles 
with  easy  axis  aligned  perpendicular  to  the  narrow  dimension  of 
film,  and  m’Lh  uniaxial  anisotropy  field  less  tlian  4;d.L  The 
magnetization  gets  -stuck*  in  a  random  .  but  essentially  unchanging 
domain  wall  configuration. 


:Z2.7 


^cdJ.  ^e^e..Y  t;, 

/  ^  cLont^in.  Iv^tL, 

(j/7ST/)5L&  f|  >  (^ 


VI.  Observations:  Amoldussen  &  Tang,  IEEE  Trans.  Mag.  22,  889,  1986. 
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Finite  Size  Scaling: 


SCALING  OF  THE  TAIL  OF  THE  DISTRIBUTION. 

THE  UNIVERSAL  SLOPE  PERSISTS  OVER  TWO  TO 
THREE  DECADES,  BUT  AT  THE  HIGH  END  (RARE, 
LONG-  LASTING  AVALANCHES)  THE  CURVES  FOR 
DIFFERENT  ANISOTROPIES  AND  SAMPLE  SIZES 

PART  COMPANY.  HOWEVER,  WE  CAN  SCALE  THE 
RESULTS  WITH  RESPECT  TO  THESE  PARAMETERS 
SO  THAT  THE  CURVES  ONCE  AGAIN  COINCIDE 
/approximately.  (THIS  IS  KNOWN  AS  FINITE- 
SIZE  SCALING.)  IN  THE  MOST  FAMILIAR 
PROBLEMS,  THERE  IS  ONLY  ONE  PARAMETER,  I  HE 


IN  THE  PRESENT  CASE,  WE  HAVE  TWO- 
PARAMETER  SCALING  ,  WITH  RESPECT  TO  THE 
LENGTH.  L,  OF  TAPE  SEGMENT,  AND  S,  THE 
CRITICAL  SLOPE. 

WE  FIND  THE  SCALING  LAW 

P(N:S,L)  =  ^g2(nS“^:LS'^) 
g2(x;y)  =  x-«  e" 

AND  B  AND  C  ARE  CONSTANTS  AND  X  =  cx^, 
S=.g/e  .  THE  INDEPENDENT  CRITICAL  INDICES  ARE 


Difficulty  of  complete  DC  erasure  of  signal  using  the  same  algorithm  as  in 
construction  of  the  zigzag: 

Initial  signal  =  reversed  magnetization  between  the  two  dotted  zigzag  lines 
Final  signal  =  remaining  lozenge-shaped  island  after  application  of  very 
large  dc  field  opposing  the  magnetization  direction  between  the  two  zigzags 


correlation  pattern  forms 


High  power  ferromagnetic  resonance, 
presence  of  thermal  noise. 


Approach  to  new  fixed  point  in 


XIIL 


The  phase  of  each  excited  spin  wave  b  is  random, 
except  that  it  is  tied  to  the  phase  of  b,  by 
the  phase  of  the  signal. 

The  spatial  autocorrelation  function  of  the  transverse 
magnetization  along  the  direction  (z)  of  the  field 
is 

c(z)=/m+(f +ez)m“(r)j=2&^ 
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Forward  Volume  Wave  Envelope  Solitons 
in  a  Thin  Unpinned  Y  I  G  Film 


Mincho  A.  Tsankov,  Ming  Chen,  and  Carl  E.  Patton, 

Department  of  Physics,  Colorado  State  University,  Fort  Collins,  CO  80523 


Workshop  on  Nonlinear  Interaction  in  Magnetic  and  Magnetooptic  Materials, 
December  12-14,  1993,  Costa  Mesa,  CA 
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Earlier  works  on  FVW  solitons: 

Kalinikos  a/. ,  Pis’ma  Zh.  Eksp.Teor.  Fiz.  38,  343  (1983) 
(pinned  film;  first  observation  of  MSW  soliton  in  YIG) 
Kalinikos  et  aL,  ffiEE  Trans.  Magn.  26,  1477  (1990). 
Kalinikos  et  aL,  Phys.  Rev.  B  42,  8658  (1990). 

(unpinned  film;  observation  of  single  and  multiple  solitons) 
De  Gasperis  et  aL,  Phys.  Rev.  Lett.  59,  481  (1987). 

De  Gasperis  et  aL,  J.  Appl.  Phys.  63,  4136  (1988). 
(investigation  of  peak  pulse  powers) 


The  Present  Work: 

-  Verifies  and  comments  on  earlier  experimental  results. 

-  Provides  some  new  results: 

-  investigation  of  the  pulse  decay  rate  vs.  power, 

-  excitation  of  solitons  in  wide  frequency  band, 

-  soliton  collision. 


_ 


/  / 


GGG 


YIG, 


GROUND  PLANE 


DIELECTRIC 


Z.AX. 


TRANSMISSION  LOSS  (dB) 


The  Experimental  Setup 


Transmission-loss  Curve:  Dispersion  Curve. 


FREQUENCY  (GHz) 


OUTPUT  PULSE  POWER  (mW) 


Second  Experiment  (A) 

Variable  Pulse  Width,  Constant  Pulse  Power 


INPUT  (Pin  =  2  W) 


OUTPUT 


Third  Experiment 

Excitation  in  a  Wide  Frequency  Band 


T  (7(3=  12  ns) 

✓ 

f  =  5.8  GHz 


OUTPUT 


n  1-5 
1.0 
0.5 
0.0 
1.5 


I  t  I 


f  =  5.7  GHz 
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Fourth  Experiment  (A) 
Pulse  Decay  Rate  vs.  Power 


TIME  (ns) 


(mW) 


Fourth  Experiment  (B) 
Pulse  Decay  Rate  vs.  Power 


INPUT  PULSE  WIDTH  (ns) 


INPUT  PULSE  POWER  (W) 


Fifth  Experiment 
Soliton  Collision 
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Theoretical  Background 


The  Nonlinear  Schrddinger  Equation: 


+V— +  77M)+-a);^-iV|u|"M=0  , 

^  dt  ^dx  2  dx^ 


_  m{x,t) .  din 

dk 


dv^  ^ 

"  g  ^ 


^k  = 


CO 


dk  dk^ 


N  = 


do) 


u 


The  necessary  condition  for  “bright”  soliton  solution: 

<  0  ;  N>0 


6) -  <<>„=  r(H^^ -4nM^)  =  - 1  h  p)] ; 


awp 


=  r4'rtMj=o„ 


2^2, 


The  analytical  solution  of  NSE: 


u(x,t)  sech 


WoC 


-2rjt 


N 


{kx-  Qt) 


Condition  for  multisoliton  excitation: 


7;>(n-V2) 


TU 

/  / 

v.ikH 

(n  -  V2)  •  C 
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Discussion 


Conclusions: 


Old  experiments  were  verified  and  extended. 

Soliton  attenuation  rate  was  confirmed  to  be  close  to 
Solitons  were  excited  in  a  wide  frequency  band. 
Collision  of  solitons  was  observed. 


E?XiSTENCE  OF  SPIN  WAVE  SOLiTOMS  !N  AN 
ANTIFERROMAGNETiC  FILM 


Jj^  * 

A  D  Boardmanrs  A  Nikitov  and  N  A  Waby. 


*j.ouie  Laboratory,  . 
Department  of  Physics, 
University  of  Salford, 
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•  :l;Rtir0d!UGtio:n 

•  B'riigiht  soliton  conditions 

•  Sbiirtace  waves  on  ain  antiferrpmagnetic  film 

isGroup  velocity  dispersion 
»  Nonlinear  coefficients 

•  ^otairriie  waves 

»  Tangentially  magnetized  ^ 

■  Perpendicularly  magn 

•  Conclusions 


o  magnetostatic  envelope  soiitons  in  thin 
ferromagnetic  films  have  received  a  significant 
deg;ree  of  attention, 

®  higih  guality  antiferromagnetic  films  made 
by  mo'lecular-beam^epitaxy  are  now 

®.antiferrom  films  may  have  many  new 

applii  cations. 

®  antiferromagnetic  resonance  frequencies  and 
tbe  frequencies  of  excited  spin  waves  belong 
to  the  infra-red. 

®  this  fact  makes  the  use  of  antiferromagnetic 
media  for  applications  very  attrqctiye. 

®  here  we  discuss  the  formation  of  spin  wave 
envelope  soiitons  for  different  combinations 
of  external  magnetic  field  and  spin  wave 
propagation  directions. 

®  Reference:  A  D  Boardman,  S  A  Nikitov  and  N  A 
Waby.  Phys.  Rev.  B48,  1  November  1993. 


■•Use  coupled  set  of  ttae  IVlaxweH  and 
Lan  dau-Liftsh  itz  eqiU;atiions . 

“Motnfeear  Schrodinger  Equation 

'M  =  isgn(R  )3^  -'Y  sgnl^f)  |  A  f  A 
dz  2  " 

®  A-:  ai^iP'Ktude  of  the  ^  ^ 

ern^/etepe  of  the  pulse  I  /1i>  t 

®T;  tfee  measured  in  a  frame  of  reference  that 
moves  with  the  pulse 

/?\  1^  -  ^  ui  ^  %  #  .1^  I  1  +\  /  ^  I  KOI  lO 


®  (B  = 

U2 


2  igroup-velocity  dispersion 


•  SiginifS.)  :sign  of  B, 


bk 

2(1  A 


:nonlinear  coefficient 


sgnCY)  :sign  of 

Necessary  condition  for  the  formation  of 
bright  solitons  is 


6,'Y<0 


Li ghthiil  Criterion) 


jSunace  Waves  on  an  Antiferromagnetic  Him| 


The  dispe.r*sJ,Oii|)  reJaitiiOn 


exp.(2;kd) 

lu^ami  jUjare  the  magnetic  permeability 
lisilisor  components  :  v  ?y 


l^v= 


'2  •  2  2  • 

1  +  R  g  (oo+oo_-  (jo  ;  j 

2  2\ /  n  ZA-txZv 

-00  )(00  -  00  ) 


_  RV(cO--Cj0+)00 

-  -  (w^oo')(co'-oo"l^Y  , 

:  t'f:  i'R^^sT m,Hy  : '1  . 

■  X  :exchange  constant 

■  M  :NET  magnetization  ®  B.B-,  •^'^'sotropy 

constants 

■  g  igyromagnetic  ratio 


iSurface  Wavesj 

®  If  Ho  =  P  till#  ,S;Unfa,ce  p-ropagate  Wfiitlnltv] 

t'fna  t^j©  fjneigiUiepiey  r©ig;lo!nrS 

.  gHc<  w  <g(H'.+R")’'' 

2 

O  2  0  0  ^/2 

g'(Hc.+^ )  <  jA>  5.9  C 

veloGity  dispersion 

B,=d'k 

c)C0^ 


:ic>3 


I  Non!inear 
[ . 


*  ys.e  method  of  A  K  Zvezdin  and  A  F  Popkov 
[iS.ov.  Phiys.  JETP  350  (1983)} 

»  fpr  smai  ac  deviations  (due  to  the  wave) 
feim  equilibrium 


•  »  the  limit  kd<<1  j  A  | 

ffl  A  is  dimensioniess  amplitude  of 
magnetostatic  potential 


®  niOitninear  coefficients  for  the  lower  and 
upper  branches  are 


*  Lighthill  criterion  satisfied  for  upper  branch 

®  Lighthill  criterion  not  satisfied  for  the  lower 
branch 


Tange ntiaily  Magnetized  Antiferromagnetic  Film  | 

1 

*  diispersion  relation 

tanJkd//|U^;)-,:^^,^ 

®  spn  wave  spectrum  (for  Hoy=0)  consists  of 
two  branches 

g(Hj,+Ho)<C0  <C0^  Lower  Branch 
9(Hc-Ho)  <  (tJ  <0)2 


2  •  2  2  2  2 


,  4_4  4  2  4.  2..  2y2,y2 

(g  R  +8R  g  H^+16g  H^H^)  ) 


ZC8 


J 


iVoSume  Wave  si 


Tan g e nti al! y  Svl agn e t i zed  Antiferromagnetic  nlm 
,  Ho=0:  only  one  branch  exists  with 

anomalous  dispersion  frequency  region 

gH^<co<g(H2+R  ) 

•  ©rbup  Velocity  Dispersion  ■ 

2^2  ,Z 

5,.i  aRcL^<o 

Vt  4g(H  +R')'‘ 


®  MiroiUiD 


Velocity 


V.=  -g'n'd 

2(g^H'+R^f 


fnear  Coefficient 


Y=— 

▼  \  f 


2 

^2  (2H,H  +H")’ 


criterion  fulfilled 


Z70  • 


IVolume  Waves! 

Perpendiculariy  Magnetized  Antiferromagneticj 
Film  _  _ _ _ _ 

•  GiiiSip.eirsi.on  relation 

Symmetric  modes  cot  /ju^kd 

Antisymmetric  modes  tan  V/J^  kd  . 

•spectrum,  again,  consists  of  two  branches 

Lower  branch  (jO^<CO  <g(H^-H^) 

Upper  branch  g(Ho+Ho)<00<Cp2  . 

•  zero’  external  field  (Hq.=0):  only  pb®  branch  in 

gH<,<  CO  <g(Hc+R^)''' 


-..,.-<-'..‘.r.o,  — . ^  . . -  - 

• 

Volume  Vv'aves 

I 

i 

] 

• 

Perpendicularly  Magnetized  AntiferromagnetiQ 

Film  i 

•  ( 
A 

pOiUip-veloci'ty  dispersion  for  this  (Hq=0) 

9  -| 

dranich  Is 

• 

1  gH,d" 

GVD 

.  tGroup  .  ’ , 

• 

.  :..u2'  ■  ■ 

,  Velocity  . 

•  1 

• 

nioniinear  coefficient 

Y  _  2  gHeHa^j-, 

V  V3  ^ 

• 

o 

LlgihthlH  cr iterid  n  n  ot  sat i  sf ied 

■  A  ^  '  ’  .  V 

Conclusions 


•  Envelope  solitons  of  dipole  spin-wave  type  are 
possible  in  two-suiblattice  uniaxial 
antiferromagnetic  films. 

®  If  the  external  magnetic  field  is  parallel  to  the 
axis  of  anisotropy,  and  is  in  the  plane  of  the 
film,  surface  and  volume  dipole  spin  waves 
«pro:pagate  within  the  film. 

®  ;lf  the  external  magnetic  field  is  perpendicular  to 
the  film  surface,  and  parallel  to  the  anisotropy 
axis,  volume  dipole  spin  waves  can  exist  in  the 
antiferromagnetic  film,  for  rather  weak  fields  (H^ 
<Hc). 

•  For  the  tangentially  magnetized 
antiferromagnetic  film,  the  necessary  and 
sufficient  conditions  for  envelope  soliton 
form 

•  Hence,  solitons  can  exist  for  the  surface 
waves  belonging  to  the  upper  branch  of  the 
spectrum.  Solitons  can  exist  for  volume 
waves  belonging  to  both  branches  of  the 
Spectrum. 

•  It  is  explained,  simply,  why  solitons  cannot 
exist  for  volume  waves  propagating  ih  a 
normally  magnetized  antiferromagnetic  film. 
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Controlling  Chaos  in  Nonlinear  Magnetic  Films 

Outline 


I  Chaos  in  Bulk  Ferromagnets 

II  Chaos  in  Thin  Film  YIG  Disks 

III  Controlling  Chaos 

IV  Concluding  Remarks 


Disigram  of  the  experiineiit<d  setup  used  for  nonliuear  ferro¬ 
magnetic  resonance.  The  magnetic  field  is  supplied  by  a  Varian  15 
magnet  controlled  by  the  magnet  controller.  X-band  microwave  power 
is  supplied  by  a  Varian  X-13  klystron  which  produces  about  50  mW  of 
power  at  about  9.2  GHz.  After  passing  through  an  isolator,  precision 
attenuator,  circulator  and  slide  screw  tuner,  the  microwaves  are  inci¬ 
dent  on  a  TEio2  cavity.  The  cavity  is  nearly  critically  coupled  so  only 
a  small  amount  of  microwave  power  is  reflected  from  the  cavity  at  its 
resonant  frequency.  The  reflected  power  is  routed  via  the  circulator  to 
a  crystal  detector,  and  from  there,  the  absorption  signal  is  analyzed  ^<gO 

electronically. 


Perpendicular  Resonance 
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Ferromagnetic  resonance  spectra  at  varying  microwave  power.  The  microwave  power 
increases  from  the  bottom  spectrum  to  the  top.  The  fuzzy  bulge  at  the  high  fre¬ 
quency  side  of  the  spectra  indicate  an  instability  in  the  microwave  absorption  called 
auto-oscillation. 


Bifurcation  Route  to  Chaos  in  the  FMR 
Absorbtion  Signal  as  the  Driving  Field  h 
is  Increased 


Time  (microseconds) 

As  the  amplitude  of  the  microwave  driving  field  h  is 
increased,  the  period  doubles  (bifurcates),  and  at  some 
critical  amplitude  of  the  driving  field  (h^),  the  sample 
absorption  becomes  chaotic. 


Magnetic  materials  display  a  rich 
variety  of  nonlinear  dynamic  phenom¬ 
ena. 


•  Regular  auto-oscillations 

•  Period  doubling 

•  Irregular  period  oscillations 

•  Intermittency 

•  Chaos 

•  Periodic  windows 


Investigations  centered  on  magnetic 
systems  are  motivated  by  several 
advantages: 

®  Magnetic  resonance  is  a  well  established  field. 

0  The  chaotic  response  occurs  in  a  homogeneous 
bulk  media. 

#  At  resonance  the  precession  angle  of  the 

magnetization  vector  increases  with  a  driving 
parameter. 

#  Nonlinear  effects  generate  both  temporal  and 

spatial  fluctuations  in  the  media. 

#  The  resonance  occurs  at  GHz  frequencies. 

0  The  auto  oscillations  occur  at  0.5-1 5MHz. 

#  A  number  of  parameters  can  be  readily 

controlled: 

magnetic  field 
frequency 
temperature 
geometry 

material  parameters 


# 


0 


0 
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Spin  0-£aniiItonian 


-  K^isX  +  SX  +  SX) 


D^SX  -is.r.Syr^) 


2J^S, 


2 

'■.y 


=  -(i>  +  b')i2S-b'byi^ 
Sy  =  '  (b-b')(2S-b'-by^ 

jU 

=  S-b'b 


.The  Hamiltonian  involves  terms  associated  with  the  interaction 
of  the  sample  with  the  static  field,  the  anisotropy  enery,  the  dipole- 
dipole  interaction  energy,  and  the  exchange  interaction  energy.  The 
spin  operators  are  expressed  in  terms  of  Bose  raising  and  lowering  op¬ 
erators. 


‘H 

•  •  •  • ) 

by.  +  ±)jJ  satisfy  3^arailton''s  Equations 


db 

dt 


ai); 


The  low  power  spectra  is  equivalent  to  neglecting 
all  terms  beyond  those  that  are  quadratic  in 


This  determines  the  normal  modes  of  the  system 


Higher  order  terms  in  bj^  add  nonlinearities  and 
chaos 


XSG 


In  terms  of  bose  operators 


=  j:<^khkh 


* 


and  the  spinwave  dispersion  relation  is 


~  ujmN,  +  uj^{akf  +  ua  + 


1 2 


2  sin^(^) 
4 


•  ^2  =  JHq 

•  =  47tjM 

•  UJ A  =  2aKA/M  where  a  is  the  con¬ 
tribution  of  the  anisotropy  energy 
depending  on  the  crystaUograpHc 
direction  of  the  applied  field,  Hq. 


The  frequency  values  above  are  for  a  ferromagnetic  sphere.  The  uniform  resonance 
of  the  sphere  ^  Woximately  10'“  spin  wave  modes 


2-S7 


Nonlinear  terms: 


Three  magnon  terms: 

=  li:Aokbihb.k. 

^  k 

Allows  the  uniform  precession,  6o,  to  exite  ±k  Suhl 
spinwave  pairs  in  the  2nd  order  Suhl  instability. 

Four  magnon  terms: 

In  the  random  phase  approximation,  the  exited  ±k 
spinwaves  are  coupled  to  degenerate  spinwaves. 

E  Skybkb^kbl,b*_y  +  c.c. 

The  terms  renormalize  the  spinwave  frequei 
cies, 

=  <^k  2J2  Tk  k'  I  by  P  . 

k' 

^k,y  terms  anmhilate  pairs  of  db/s  spinwaves  while 
creating  ±k'  pairs. 


=  E  T,,yb,blbi,,bl,+l 

k,k'  Z 


In  a  sphere,  there  are  10^®  degenerate  spinwaves! 


2.'^8 
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Truncating  the  number  of  active  modes  still  yields  qualitatively  good  results. 
On  the  left  are  signals  from  a  truncated  model  compared  to  experimental 
signals  on  the  right.  From  M.  Warden,  thesis.  Despite  the  large  number 
of  degenerate  modes,  experimental  chaotic  signals  often  have  fractal  dimen¬ 
sions  between  3  and  5,  indicating  that  the  system  reduces  the  number  of 
degrees  of  freedom.  This  suggests  that  cooperative  modes  are  developed. 
Intensive  effort  is  underway  to  determine  the  nature  of  these  normal  modes 
in  bulk  magnetic  materials. 


-4  -2  0  2  4 
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Unlike  FMR  in  spheres  where  the  resonance  is  degenerate  with  a  large 
number  of  spin  wave  modes  in  the  middle  of  the  spin  wave  band,  in  a 

disk,  the  driven  modes  are  at  the  bottom  of  the  band,  and  relatively  few  9 

modes  are  excited.  This  is  due  to  the  difference  in  the  demagnetization 
factors  between  a  sphere  and  a  disk. 
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For  a  circular  film,  the  spatial  functions  are  of  the  form 


P I 


sin(z/(;^) 

cos[u(p) 


where  x,,,,  is  a  zero  of  the  Bessel  function  and  a  is  the  sample  radi 


Direct  Modes 


Hidden  Modes 


Domains  of  the  transverse  component  of  magnetization  of  the 
lowest  magnetostatic  modes  in  a  circular  disk. 
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Figure  l  FMR  absorption  and  the  domain  of  oscillatory  behavior  over  three 
decades  of  pumping  power.  Oscillatory  behavior  is  indicated  by  horizontal 
stripes  at  IdB  intervals  in  the  real 
oscillation  arises  from  each  mocTe 
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The  Hamiltonian  for  a  disk  includes  the  interaction  with  the  static  field,  the  de¬ 
magnetization  field,  the  dipole-dipole  interaction  and  the  microwave  pumping  field. 

The  magnetization  operators  are  expanded  in  terms  of  the  normal  modes  of  the  disk 
(Bessel  functions).  From  the  Hamiltonian,  the  equations  of  motion  of  the  mode  am¬ 
plitudes  can  be  derived.  Note  that  the  non-linear  term  arises  from  the  dipole-dipole  2.94^ 


By  including  the  (2,1)  and  (1,1)  hidden  inodes  along  with  the  (0,1), 
(0,2)  and  (0,3)  direct  modes,  “thumbs”  seen  in  the  experiment  can  be 
reproduced,  suggesting  that  the  hidden  modes  play  a  role  in  the  system 
dynamics. 


frequencies  (MHz) 

Onset  power  (dB) 

Finger  # 

model 

expt. 

model 

expt 

1 

3.97 

4.61 

-14.5 

-13.5 

1.5 

1.5 

-13.7 

-12.5 

2 

— 

5.22 

— 

-12.5 

7.8 

9.2 

-14.0 

-12.5 

Table  5.2:  Auto-oscillation  frequencies  CcJculated  from  the  five  mode  model  and 
values  obtained  from  sample  1. 
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CONTROLLING  CHAOS 


# 


Goal  -  stabilize  certain  unstable  orbits 

Three  recognized  techniques: 

1.  Ott,  Grebogi,  Yorke  (OGY)  PRL  ’88 

Discrete  method  of  applying  a  perturbation  at 
intervals  determined  by  behavior  in  phase 
space. 

2.  Hunt  PRL  ’91 

Discrete  method  of  applying  a  perturbation  at 
intervals  determined  by  behavior  in  time 
space. 

3.  Pyragas  PL  A  ’92 

Analog  method  applying  a  continuous 
"feedback"  perturbation. 
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OGY  Method  -  discrete 


Consider  a  path  in  N-dimensional  phase 
space 


Goal  is  to  stabilize  the  orbit  such  that  it 
regularly  penetrates  the  plane  at  a  fixed 
point  Op 


A  region  for  O^  near  Op  is  defined  over 
which  a  linear  response  to  a  perturbation 
of  some  systematic  variable  is  expected 


300 


Behavior  is  expanded  in  a  region  near  Op  for 
which  a  linear  approximation,  valid  for 


On  -  Op  <  AO 


then 

^n+m  “^F  ^ 


M  is  a  2x2  matrix  with  eigen  vectors  along 
the  stable  manifold  and  the  unstable 
manifold  (saddle  point)  associated  with 
the  perturbation 


Figure  3.  Outline  of  OGY  control  algorithm  for  a  saddle  fixed  point:  (a)  the  n^i 
iterate  falls  near  the  fixed  point  ^f(po)-  (h)  Turn  on  the  perturbation  of  p  to 
move  the  fixed  point,  (c)  The  next  iterate  is  forced  onto  the  stable  manifold  of 

5f(po)-  Turn  off  the  perturbation. 


From  Ditto,  Rauseo,  and  Spano. 


Controlling  Chaos  in  a 
Magnetoelastic  Ribbon 


The  experimental  setup  of  Ditto,  Rauseo  and  Spano  who  used  the  OGY  method  to 
stabilize  the  chaotic  swaying  of  a  magnetoelastic  ribbon.  The  ribbon’s  Young’s  mod¬ 
ulus  was  highly  sensitive  to  external  magnetic  fields,  so  the  external  AC  magnetic 
field  induced  chaotic  motion.  Ditto,  Rauseo  and  Spano  were  able  to  apply  discrete 
perturbations  to  the  magnetic  field  to  create  a  periodic  motion  in  the  ribbon  under 
circumstances  that  would  be  chaotic  without  perturbation. 
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1000  2000 

Iteration  Number 


m 


The  ability  of  DRS  to  svsntch  the  periodicity  of  the  ribbon.  Notice  the  brief 
chaotic  transients  between  switchings.  From  the  Proceedings  of  the  1st  Experimental  Chaos 
Conference. 
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2-®  2.5  3.0  3.5  4.0  4.5  5.0 

Previous  Position 


Fig.  4.  The  current  ribbon  position  Vj  is  plotted  versus  the  previous  position 
Vj.i  yielding  a  phase  portrait  of  the  attractor,  (a)  Period  one  control  (black), 
(b)  Period  two  control  (black).  In  both  cases  the  grey  represents  the 
uncontrolled  chaotic  behavior. 


2.  Hunt  Method 

A  modification  of  OGY. 

Consider  an  oscillatory  time  signal 


An  amplitude  window  W  establishes  the 
application  and  magnitude  of  the 
perturbation 

F  =  L[A„-W^]  5 

5  =  lif  <W/2 

0  if  Ajj  -  >  W/2 

Example  is  a  resonator  circuit  driven  at  53 
kHz 

Chaotic  signals  have  been  stabilized  in 
periodic  orbits  up  to  period  86! 
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First  return  map  (I„+i  vs  I^,  left 
By  applying  the  control  signal  i 
(upper  graphs)  and  period  21  (Ic 


time  ( ms ) 


current  through  chaotic  oscillator  (right), 
in  (c),  the  system  stabilized  on  period  5 
hs)  orbits.  3*97 


3.  Pyragas  Method 
An  analog  technique 

F(t)  is  a  perturbation  to  an  external  control 
parameter 

F(t)  =  K  [yi(t)  -  y(t)]  =  K  D(t) 

yj(t)  is  the  desired  behavior  (stable  orbit) 
of  the  system 

y(t)  is  the  actual  response 

D(t)  is  the  difference 

K  is  the  adjustable  weight/drive  factor 

Example  is  Rossler  system  and  the  Lorenz 
system 
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Figure  2:  The  output  signal  y(t)  and  the  forcing  term  F(t)  of  the  fal  Rossler  system-^ 7^/'.  = 
~y  ~  <fy/ dt  =  X  +  0.2y  +  F(t),  dzjdi  =  0.2  +  2(1  —  5.7),  stabilized  to  a  period-5  cycle  and 
the  TiOrniB  5Ystem:<fz/dt  =  10(i  - y),  dyfdi  =  —xz  +  28z  -y  +  F{t),  iz/dt  =  xy-  (8/3)2, 
stabilized  to  period-2  motion.  The  inset  is  an  x-y  portrait  of  the  stabilized  system. From 
Pyragas 


// 


Figure  3:  The  graph  above  shows  the  leAding  Lv^punov  exponent  A.  a  measure  of  the 
system’s  dependence  to  initiaj  conditions,  versus  K  in  stabilizing  orbital  periods  of  1  and  2 
of  the  Rossler  system.  When  A  is  less  than  zero,  the  system  is  no  longer  sensitive  to  initial 
conditons.  From  .Pyragas. 
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Controlling  Chaos  in  YIG  Film 
using  Feedback  Modulation 


Magnetic  field  Ho  is  modulated  by  the  chaotic  signal  S(t)  obtained 
from  the  FMR: 


Hq  — *•  Hq  +  ^H[S(t  —  t)  S{t  t)] 


AH:  perturbation  amplitude 
r;  delay  time 

S{t):  DC  component  of  FMR  signal 

The  perturbation  term  is  proportional  to  the  AC  component  of  the 
FMR  signal. 


Diagram  of  FMR  experiment  for  closed  or  open  loop  feedback.  For  open  loop  feed¬ 
back,  the  function  generator  is  engaged  to  perturb  the  sample,  while  in  closed  loop 
feedback,  the  delay  line  is  employed.  ^ 
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^  Experimental  results  of  controlling  chaos  of  FMR  using  the  time  delayed  feedback  method,  a) 

Chaotic  oscillation  of  FMR  without  perturbation  converted  to  b)  period  4  (gain=0.45  arbitrary 
unit),  c)  period  2  (gain=0.80  a.u.),  d)  period  1  (gain=0.50  a.u.),  and  e)  the  quiescent  state 
(gain=6.50  a.u.),  as  the  gain  of  the  amplifier  wa.s  increased. 
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of  the  conversion  of  a  p 
nal  to  a  cyclical  plot  in  phase  s 
y=sincjt 

- - - - 1 - - - 1 


Example  of  the  conversion  of  a  period  two 
time  signal  to  a  cyclical  plot  in  phase  spac 

y=sincjt+0.3sin[(w/2)t  +  7r/ 4] 


N  dimensinal  phase  space: 
y  =sin[cj(t+nT)  +  0.3sin[aj/2(t+nT)+7r/ 4] 

n=0,  1,  2 .  N-1 


An  example  of  a  2— D  phase  plot 
with  cjt=7t/3 


The  double  loop  cycle  represents  a 
period  two  response  (bifurcation) 


Computer  Simulation  of  Controlling 

Chaos 


The  normal  modes  of  FMR  in  circular  disk  are  the  magnetostatic  modes  with  the 
space  distributions  of  Bessel  functions.  In  the  computer  model,  the  magnetization 
of  YIG  film  is  expanded  as  the  superposition  of  these  normal  modes  with  the  time 
dependent  amplitudes,  a.-.  The  equations  of  motion  of  can  be  derived  from  the 
Hamiltonian  of  the  system 

^  =  -i7  ({Ho  -  Hr  -  *T)a,-  +  y/r  +  27rM,  ^  Ajwa'aiaij  , 


Here  is  the  low  power  resonance  field  of  the  ith  mode,  hp  is  the  rf  pumping 
The  coefficients  Aijki  introduce  the  nonlinear  coupling  between  these  modes.  e 
nonliear  dynamics  of  high  power  FMR  in  thin  cirlular  YIG  films  can  be  simulated 

by  this  computer  model. 


Time  Delay  Plot  of  an  Experimental 

Chaotic  Signal 
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Time  Delay  Plot  of  a  Chaotic 
Experimental  Signal  Controlled  to 
a  Period— 4  Signal 
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Time  Delay  Plot  of  a  Chaotic 
Experimental  Signal  Controlled  to 
a  Period-2  Signal 
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AC  Absorption(t  +  T) 

t  =  50  ns 
Gain  =  0.8 

Feedback  Delay=470  ns 
Arbitrary  Absorption,  Gain  Units 


Time  Delay  Plot  of  a  Chaotic 
Experimental  Signal  Controlled  t 
a  Period- 1  Signal 
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Time  Delay  Plot  of  the  Transient 
Response  of  the  Controlled 
Chaotic  Signal  to  the 
Quiescent  State  in  the  Model 
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System  evolves  chaotically  for  2  microseconds  before  the  feedback 
is  turned  on,  and  after  a  transient  period  of  about  5  microsends, 
the  system  is  controlled  to  the  quiescent  state.. 
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Computer  simulation  results  of  controlling  chaos  of  FMR  using  the  time  delayed  feedback 
method,  a)  chaotic  signal  without  perturbation.  With  increasing  gain  of  the  feedback  loop  Hp, 
chaos  was  converted  to  b)  period  4  osciUation,  c)  period  3,  d)  period  2,  e)  period  1,  f)  quiescent 
state.  The  perturbation  was  turned  on  at  t=0.  g)  The  system  returned  to  the  chaotic  state 
after  turning  off  the  perturbation  at  t=100ns.  3  2  S’ 


Modulation  field  [G] 


/e  cp 

AHCO  -  AH 


Modulation  frequency  [kHz] 


Inside  the  triangular  regions  in  the  modulation  field-frequency  plane,  the 
chaotic  FMR  signal  was  stabilized  to  periodic  orbit  by  modulating  the  mag¬ 
netic  field  with  a  sine-wave  perturbation.  The  numbers  denote  the  frequency 
(in  MHz)  of  the  periodic  oscillation  after  stabilization. 
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CONCLUSIONS 


Nonlinear  dynamics  in  magnetic  materials  has 
been  and  continues  to  be  an  exciting  and 
active  field. 


Numerical  models  of  the  equations  of  motion  are 
readily  available  from  the  0-£amiltonian. 


Thin  film  geometries  limit  the  number  of 
degenerate  modes  and  truncation  is  not  a 
serious  limitation. 


By  appropriate  perturbation  techniques  it  is 
possible  to  return  the  system  from  a  chaotic 
state  to  the  quiescent  state. 


Chaos  can  be  Controlled 
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